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1 . INTRODUCTION 


This General Test Plan is submitted to NASA MSFC as the second major 
deliverable on Contract NAS 8-27335. The plan defines the test series 
necessary for the completion by MSFC of a redundant sensor strapdown IMU 
evaluation experiment. 



1.1 EXPERIMENT OBJECTIVES 

This experiment series uses a prototype (breadboard) strapdown IMU. 

It contains redundant gyros (six in a dodecahedron), and three orthogonal 
accelerometers. The system software incorporates failure detection and 
correction logic and a land vehicle navigation program. 

A significant feature of this experimental program is the use of 
Single Axis Reference units (commonly called single axis platforms) for 
attitude reference. This device is described in Reference 1. Hardware 
experimentation with redundant strapdown systems using conventional 
single degree of freedom gyros has been done elsewhere (MIT, ERC). 

The principal objective is a demonstration of the practicability, 
reliability, and performance of a SAR redundant strapdown IMU system 
with failure detection and correction’, in operational environments. 

Engineering information will be obtained on permissible failure 
threshold settings with real instruments and the accumulation of 
navigation errors before a failure is detected. 

1.2 SCOPE 

This plan includes: 

• Equipment Configuration 

• Error Analysis 

• Software Requirements 
t Test Descriptions 

Configuration 

The equipment configuration a'nd software requirements were estab- 
lished in the Program Definition and Experiment Configuration Plan, 
Reference 2, and only a summary and revisions are included in this plan. 

Error Analysis 

Design of the system tests and analysis of the test results requires 
analysis and modeling of system errors. A preliminary error analysis of 
the MSFC Breadboard Dodecahedron SAR Cruise Navigation system has been 
included in this plan. The analysis is in sufficient depth to allow: 
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1) evaluation of calibration and compensation requirements 

2) test design (sequence, input parameters, duration) 

A continued error analysis effort will be required by MSFC for inter- 
pretation of test results, setting performance criteria, and refinement 
of the test sequence. 

Test Plans 

The test objectives, constraints, description and analysis 
(Section 5 through 9) are divided into five test phases, as shown in 
Figure 1-1. The test descriptions are based on the assumption that the 
hardware and software configurations described in Sections 2 and 4 have 
been implemented and completely checked out. The test sequence is 
arranged such that the calibration, alignment, and navigation functions 
are completely checked out prior to exercising the FDDC logic with 
failures, simulated failures, and degraded operation. 
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Figure 1-1. Overall Test Sequence 
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2. EQUIPMENT CONFIGURATION 


The configuration of the test article, test station, and mobile van 
upon which the- test series is based are documented in Reference 2. Only 
changes and additions are listed below. 

2.1 BREADBOARD DODECAHEDRON (BB DDH) 

The accelerometer torquer loops are being modified to reduce the 
quantization from .4 Fps/pulse to .15 Fps/pulse. This will reduce the 
alignment filter convergence time and reduce the noise in the calibration 
and navigation routines. 

The SAR which did not have a gimbal reference position output pulse 
is being reworked. All six SAR's will have a zero reference pulse. 

Modifications are being made to the BB DDH electronics package to: 

1) Combine positive and negative output pulses from each sensor 
onto a single output line to reduce the number of recording 
channels required in the van. 

2) Combine three fault monitors (wheel current, gimbal spin 
detector, and encoder error) into a single failure discrete. 

These changes are shown in the revised BB DDH block diagram. 

Figure 2-1 . 

Another change from Reference 2, Paragraph 2.3.4 is; the three 
torquing rates available from the Control Panel (Figure 2-1) are: 

0.4 degree/hour 
4.0 degree/hour 
60 degree/minute 

2.2 MOBILE VAN CONFIGURATION 

The van instrumentation is being developed approximately as des- 
cribed in Reference 2, Section 2.4 and repeated in Figure 2-2, although 
the design is not yet complete. The alternate configuration employing 
two tape recorders in the van will not be used. Redstone has the IRIG B 
time code on a radio link and a receiver may be installed in the van. 

The van will be air conditioned. Van instrumentation is discussed in 
Paragraph 9.4.2. 
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Figure 2-1. Block Diagram BB DDH, Control Panel, and IFE Interfaces (Preliminary) 











2.3 INTERFACE DEFINITION 

The equipment interfaces are as defined in Reference 2, Section 3, 
except for the change in accelerometer quantization previously discussed. 
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Figure 2-2. Mobile Van Functional Block Diagram (Conceptual) 

















3. ERROR ANALYSIS 


A preliminary analysis of navigation errors In the laboratory and van 
tests is presented In Appendix A. In this analysis, the navigation error 
equations of a nonredundant strapdown I MU are described. This descrip- 
tion is then extended to a dodecahedron with no failures. Finally, the 
performance degradation associated with failure detection, diagnosis and 
correction (FDDC) is discussed. 

A preliminary navigation error estimate for the BB DDH is also 
presented. Performance is estimated both without and with a gyro drift 
update prior to navigation. Also, performance is estimated both without 
and with a redundant sensor program (RSP) change to compensate gyro mass 
unbalance and SAR misalignment as a function of SAR gimbal angle. 

In Appendix A, attitude and navigation error sensitivities and 
navigation error estimates are derived in Customary Units. These error 
sensitivities and error estimates are presented in Figures 3-1 through 
3-9 and in Table 3-1 in both Customary Units and in the International 
System of Units. 

3.1 NONREDUNDANT STRAPDOWN IMU 

The attitude, velocity and position errors that result from the 
following error sources in a nonredundant strapdown IMU are described in 
Appendix A. 

• initial azimuth al Ignment. error 

• initial level alignment error 

• level gyro drift 

t vertical gyro drift 

• gyro misalignment 

• level accelerometer bias error 

• level accelerometer scale factor error 

• level accelerometer misalignment 
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Error correlations involving initial azimuth alignment error and level 
gyro drift, and initial level alignment error and level accelerometer bias 
error are also described. 

It is shown in Appendix A that the attitude and navigation errors 
that result from an initial azimuth alignment error are of the same form 
as the corresponding errors resulting from level gyro drift. Similarly, 
the attitude and navigation errors that result from an initial level 
alignment error are of the same form as corresponding errors resulting 
from level gyro drift, vertical gyro drift and level accelerometer bias 
error. The latter statement is also true for gyro misalignment and level 
accelerometer misalignment. Thus, with the exception of the attitude and 
navigation errors that result from level accelerometer scale factor error 
(and with the exception of certain other second order errors), all 
attitude and navigation errors are of the same form as those resulting 
from level gyro drift, vertical gyro drift and level accelerometer bias 
error. The latter errors are plotted in Figures 3-1 through 3-9. 

It is shown in Appendix A that the attitude and navigation errors 
resulting from level accelerometer scale factor error will be much less 
than corresponding errors resulting from other error sources. Hence, 
level accelerometer scale factor errors may be ignored in this analysis. 

3.2 STRAPDOWN DODECAHEDRON 

The sensor error coefficients of a strapdown dodecahedron may be 
related to the corresponding error coefficients of a nonredundant IMU. 

Assume that the gyro bias errors of a dodecahedron are statistically 
independent with zero mean and equal variance. It is shown in Appendix A 
that for no failures, 
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Similarly, for gyro spin axis mass unbalance and assuming that the 
level acceleration input components are negligible. 
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where o g = gyro spin axis mass unbalance variance 
and the x axis is vertical. 

For gyro input axis mass unbalance and assuming that the SAR gimbal 
angles as defined in Appendix A are equal to 90 degrees, 
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where Oj = gyro input axis mass unbalance variance 

For the SAR misalignments, assuming that the level angular rate input 
components are negligible and ignoring the internal misalignments. 
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where a c = SAR misalignment variance. 

Equations extending the foregoing relations to include the effects 
of SAR gimbal angle changes are provided in Appendix A. 

3.3 FAILURE DETECTION, DIAGNOSIS AND CORRECTION 

A failure of one or more sensors in a dodecahedron will cause system 
performance to be degraded. Degradation may result from one or more of 
the following: 
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e undetected errors 


• transient errors 

• errors resulting from the deletion of sensors in the estimation 
equations 

Each of the foregoing sources of error may be evaluated with the RSP 
in the Simulation Mode and with the error equations presented in 
Appendix A. 

3.4 ESTIMATED PERFORMANCE 

A rudimentary estimate of navigation accuracy in the forthcoming 
tests is presented in Appendix A. This estimate is based on the estimated 
BB DDH performance characteristics provided by MSFC and summarized 
in Table A-l . The navigation accuracy estimate is summarized in 
Table 3-1. 

It is evident from Table 3-1 that navigation accuracy will be 
improved by the addition of a gyro drift update prior to navigation and 
by an RSP change to compensate gyro mass unbalance and SAR misalignment as 
a function of SAR gimbal angle. It appears, however, that the objectives 
of the current test series can be realized without RSP changes. (If 
necessary, a work around procedure can be devised to provide a gyro 
drift update prior to navigation). Hence, it is recommended that RSP 
changes to improve navigation accuracy be deferred. 
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Table 3-1. Navigation Accuracy Estimate 


Test Conditions 

Estimated 

Accuracy 

DDH 

rotated 

during 

navigation 

Gyro 
dri ft 
update 
prior to 
navigation 

RSP change to 
compensate gyro 
mass unbalance 
and SAR misalignment 
as a function 
of SAR gimbal angle 

n.mi ./hr 

Km/ hr 

No 

No 

No 

4.6 

8.5 

No 

Yes 

No 

1.3 

2.4 

Yes 

No 

No 

9.1 


Yes 

Yes 

No 

5.1 

9.4 

Yes 

Yes 

Yes 

1.9 

3.5 
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Figure 3-2. Velocity Error versus Time Level Gyro Drift - 0.01 deg/hr 
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Figure 3-5. Velocity Error versus Time Vertical Gyro Drift - 0.01 deg/hr 
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Figure 3-7. Attitude Error versus Time Level Accelerometer Bias — 1 0 pg 
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Figure 3-8. Velocity Error versus Time Level Accelerometer Bias - 10 Mg 
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4. SOFTWARE CONFIGURATION AND CHECKOUT 


4.1 PROGRAM MODIFICATIONS 

Reference 2 defined the Redundant Sensor Program (RSP) revisions 
resulting from the changes in configuration between ERC and MSFC 
experiments. These are differences in: 

• dodecahedron configuration 

• test equipment configuration 

• test site location 

In this section all mandatory changes are tabulated and optional changes 
are identified. Section 3, Error Analysis, establishes the approximate 
system error attributable to the optional changes. Recommendations are 
made for changes to be incorporated, and the test plans are based on 
these recommendations. 

4.1.1 Calibration Program 

4. 1.1.1 Instrument Change Discussion 

Reference 2 
Section No . Page 


Nominal Gyro Bias 

4.2.2 

4-3 

Nominal Accel erometer Bias 

4.2.2 

4-3 

Gyro Scale Factor 

No change 


Accelerometer Scale Factor 

4.2.2 

4-4 

Site Dependent Constants 

Latitude \ 

4.2.2 

4-4 

Azimuth of Assembly Axes e 

4.2.2 

4-4 

Static Test Interval Change 

Note 1 



(from 100 seconds to 240 seconds) 


4.1.2 Alignment 

4. 1.2.1 Instrument Change Section No . Page 

SIGA (quantization noise) 4.2.3 4-5 

SIG1,2,3 (initial misalignment estimate) 4.2.3 4-5 
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4. 1.2. 2 Site Dependent 


Section No . Page 


AZHT (initial azimuth estimate) 

4.2.3 

.4-6 

RILX (inertial to VEN transformation) 

Note 2 


4. 1.2. 3 Real Time Alignment Error Output 

4.1.3 Navigation and FDDC 
4. 1.3.1 Instrument Change 

Note 3 


FDDC 

• Bypass Accelerometer C matrix 

4.2.2 

4-5 

• Bypass Accelerometer FDDC 

4.2.2 

4-5 

• Accelerometer State Flag (Sa) 

4.2.2 

4-5 

• Accelerometer Nominal A matrix 

See Note 4 


• Gyro Nominal A matrix 

See Note 4 


• Gyro Nominal C matrix 

• Addition of Failure Discrete 

See Note 4 


Input Processing 

See Note 5 


Navigation 

• Sensor quantization change 

4.2.4 

4-8 

• Sensor nominal bias 

4.2.4 

4-8 

4. 1.3. 2 Site Dependent Constants 

Latitude, Longitude, Gravity 

4.2.4 

4-7 

4.1.4 Driver Routine 

• Sensor quantization, QA, QG 

4.2.4 

4-7,8 

• Sensor nominal bias 

4.2.4 

4-8 

4.1.5 Display 

Nixie Display to HP-2116 CRT Display 

4.2.6 



Note 1 : 

To change the nominal static test calibration interval from 
100 seconds to 240 seconds, change the input parameter STLT, 
static test sampling time in milliseconds to 240,000 (B23). 
(See Page 4-4, Reference 2). 

Note 2: 

Computation of the new Rj matrix is explained on Page A-13, 
Volume 1, Reference 3. 
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Note 3: 


Alignment 

If the parameters computed in Figure B63 of Volume 2, Reference 4 
are desired for a real time alignment mode, the following changes 
should be made. 

1) After the sequence number 4978 (JPL A13), insert the 
following instructions: 


JMP PD6 

2) Change the following instructions: 


Se quence No . 


Mas Change To 


4981 JMP PD6 JMP D00 

3919 JMP D00 JMP All +5 


The "K" matrix (3x3) must be initialized to the proper value at 
the beginning of a run. 

Note 4: 


1) Accelerometer A Matrix 


Z 

X 

Y 


X Y Z 


0 0 1 
0 0 0 
1 0 0 
0 0 0 
0 1 0 


2 ) 


X Y Z 


boa 
b o -a 
a -bo 
a b o 
o a b 
o a -b 


b = .52573111 
a = .8506508 
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Note 4: (continued) 

3) Gyro C Matrix (15 x 6) 


0 

0 

-a 

a 

-b 

-b 

0 

-a 

0 

b 

-a 

b 

0 

-a 

b 

0 

-b 

a 

0 

b 

-a 

b 

0 

-a 

0 

-b 

-b 

a 

-a 

0 

a 

0 

0 

-b 

-b 

a 

a 

0 

-b 

0 

-a 

b 

-b 

0 

-b 

a 

0 

-a 

b 

0 

-a 

b 

-a 

0 

-b 

b 

0 

0 

a 

-a 

a 

b 

0 

-a 

0 

b 

-b 

-a 

0 

a 

-b 

0 

-b 

-a 

a 

0 

0 

b 

a 

b 

-a 

0 

-b 

0 

a 

a 

-b 

-b 

0 

0 


Note 5: 

Failure Discrete Input Processing 

In order for FDDC to do real time discrete monitoring of the six 
gyros, some minor program changes will be required. These 
changes and the coding that already exists in FDDC and the real 
time input routine are based on the following assumptions: 

a) An 18-word list will be read in every 40MS from the IFE. 

b) The first six words correspond to gyros one through six. 

c) The accelerometer pulses are in words seven, nine, and eleven. 

d) The table angle, scaled at +15, is in word 13. 

e) Gyro monitor discrete is in word 14. 

f) The format of word 14 is as shown below: 

BIT POSITION 


012345 6789 10 11 

12 13 14 

15 16 17 

18 19 20 

21 22 22 

not used Gyro 6 Gyro 5 

Gyro 4 

Gyro 4 

Gyro 2 

Gyro 1 


GYRO DISCRETE INPUTS 
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Note 5: (continued) 

The low order 18 bits are used to store the gyro discrete data. 
Three bits are used for each gyro and if any one or more of the 
bits are on (equal to 1 ) FDDC will treat that gyro as having 
failed until all three bits are ZERO. 

If the following modifications are implemented, then word 14 will 
be saved and used by FDDC for gyro internal discrete monitoring. 

Replace the instruction JMP FO at sequence number 6596 with the 
following: 

LDA BUF1 +13 Gyro Discrete word 

STA FFDG FDDC Monitor word for gyros 

JMP FD Check for failure processing, 

begin cycle calculations 

4.1.6 Optional Program Modifications 

The program modifications discussed below appear to be desirable 
for follow-on or more advanced testing. For the initial BB DDH tests, 
defined in this test plan, it is presumed that these changes will not be 
implemented. The purposes of the BB DDH testing can be accomplished 

i 

with the current program modified as discussed in this section, and the 
procedures discussed in Sections 5 through 9. 

4. 1.6.1 SAR Gimbal Position Compensations 

The SAR instrument error model differs from the strapdown model 
incorporated in the RSP in three ways: 

1) The input axis mass unbalance drift error is a function of the 
SAR gimbal angle. 

2) There is an error term (drift) M gp , which is a function of 

acceleration along the output axis. This is also a function 
of gimbal angle. 

3) The gyro input axis misalignment is a function of the 
gimbal angle. 

The gimbal angle for each SAR is equal to the total (body fixed) input 
angle about that instrument input axis. 
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Compensation for the g sensitive drift terms would require: 

1) Accumulation of total body fixed input angle (gimbal, angle) 

p 

for each SAR, ©j. 

2) Modification of the gyro mass unbalance compensation matrix 
equation from the present form 

el = (K, U sA + 1C, U IA ) ' ^ 


- (K, a,* * ic, u IA ♦ k 3 u^) • h 

= sensed acceleration vector expressed in 
body fixed coordinates 

= gyro mass unbalance compensation parameters 

= output axis acceleration sensitivity 
= the spin axis unit vector expressed in body fixed 

p 

coordinates, a function of ©j. 

= the input axis unit vector expressed in body fixed 
coordinates, a constant vector. 

= the output axis unit vector expressed in 
body fixed coordinates, a function of ©j. 

Compensations for input axis alignment as a function of gimbal 

p 

angle would employ the same computed angles ©" in a modification of 
the equations for computing the angular rotation vector Aa, and the 
gyro test signals (V' matrix). 

The equations and coding for these changes have not been developed 
but it is presumed that suitable forms can be found that can be 
incorporated in the present DDP-124 program. 


to the form 


el 


where : a^ 


h 


“sA 


— IA 
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4. 1.6. 2 Gyro Bias Estimator 

If considerable time has elapsed since a complete gyro calibration, 
navigation performance can be improved by an update of the gyro drift 
estimate immediately prior to a navigation run. In BB DDH testing this 
estimate will be made after the fact from taped data processing using 
special procedures. 

The RSP program could be modified to estimate the total drift 
(g sensitive plus fixed drift) about each axis, in the van, using an 
optically determined azimuth. The bias estimator would contain a filter 
for gyro and accelerometer quantization noise. The filter would be 
significantly different from the self-contained alignment filter presently 
in the RSP program. The self-contained alignment filter does not 
estimate gyro drift, and it is not designed to make use of optical 
measurements . 

The size of the program required for bias estimation is comparable 
to the size of the present alignment routine. 

4. 1.6. 3 Maximum Likelihood Failure Detection and Isolation 

Recent developments in failure detection methods have shown optimal 
(maximum likelihood) filtering techniques to yield more efficient 
(reliable, accurate) system operation in the presence of degraded 
instrument performance (soft failures). 

Using these techniques, and logic matched to the model and require- 
ments of the MSFC SAR Dodecahedron navigation system, the RSP program 
could be modified for evaluation of more sophisticated and efficient 
failure isolation techniques. 

A quick analysis of a set of maximum likelihood filtering equations 
indicates that this scheme can be incorporated in place of the present 
FDDC routine. The execution time for this scheme will be approximately 
equal to the FDDC routine. The area which needs to be investigated 
further is the quantization effects of the DDP-124 fixed point 
arithmetic on some of the filter variables. 
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4. 1.6. 4 Calibration Program Changes 


If SAR gimbal angle compensations are employed, changes will be 
required in the RSP calibration routine and calibration procedures to 
measure output axis g sensitive drifts and SAR internal misalignment 
(two components) for each instrument. 

The test control logic to select rate test, static test, compute, 
etc., could be revised to make use of the revised calibration sequence 
much more straightforward (Subparagraph 5.4.2.2). 

4.2 SOFTWARE CHECKOUT 

4.2.1 Copying and Tape Manipulator 

Before starting checkout of the Redundant Sensor Program (RSP), the 
current symbolic tape and object tape should be copied and filed away. 
After a good copy has been made of the RSP tapes, the modifications 
(Subsection 4.1 above) to the program should be made and an assembly 
with a listing should be done. Once an error free assembly is obtained 
and all the modifications have been inserted in the proper place, an 
object tape should then be made. In order to do this, the FORTRAN 
Routines (RDWT and KFTR) will have to be put onto paper tape. 

In order to create a new object tape, the system DAP II Assembly 
program must be loaded into memory and the following options selected: 

Source input from tape 

Object output on tape 

Avoid rewind and End-of-file (EOF) on object tape 

Bypass listing 

The RSP tape should then be mounted on the input tape unit and a 
scratch tape on the output unit. The scratch tape will contain the 
object code. At the completion of the assembly, the RSP tape is removed, 
the system tape is mounted again and the FORTRAN compiler is selected. 

The following options should be selected: 

Source input from paper tape 

Object output on magnetic tape 
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Write an EOF and rewind object 

Tape at completion 

Bypass listing 

At the completion of the compilation, the object tape will contain the 
complete RSP except for the system subroutines. In order to load the 
program into memory, put the system tape on one unit and the object tape 
on the other. 

Select the loader with the following options: 

Magnetic tape input 

Select file one 

Put loader map on the typewriter 

At the completion of loading the program, it is possible that a system 
subroutine will be missing. This routine is 02LD anc j it is a line printer 
dump routine. Missing this routine will not cause any problems as long 
as there is not a line printer attached to the DDP-124. 

As part of the software checkout, it will doubtless, be necessary to 
make octal dumps of memory. The RSP has an octal dump routine built 
into it. In order to use it, the flag TYFL needs to be set to 1, either 
with the input routine or the assembler. Also, the octal location of the 
dump needs to be noted. 

To dump portions of memory onto the typewriter, set the location 
counter EQUAL to the octal location that corresponds to the symbolic 
location DUMP. The A register is set equal to the desired starting and 
stopping memory locations. Since the A register can only hold either 
octal digits, only the most significant address digits can be specified. 
For example, if an octal dump of memory from octal location 07010 through 
octal location 15672 is desired, first, transfer control to the DUMP 
routine, at this point the computer will halt. Next, the A register is 
set equal to the following: 

07011570, and then the start button is pushed. Assuming that the 
dump flag has been set to dump onto the typewriter (TYFL equal to 1), the 
octal dump will then take place. At the conclusion of the dump, the 
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routine returns to DUMP and halts. At this point, another dump can be 
initiated by setting the A register with another pair of addresses. 

The above A register setting will actually dump from 07010 through 
15700. Since the least significant octal number of the address cannot 
be specified, an address of 15700 was selected so that 15672 would be 
included in the dump. 

4.2.2 Checkout of Alignment and Navigation Program 

In order to do a detailed checkout of the program, a minimum of 
two runs are required. Both are nominal runs, one while in the alignment 
mode and the other in the normal mode. Both runs should be made on the 
original unmodified program and then repeated after all changes have 
been incorporated into the program . Also, the Driver/Navigation mode run 
should be made to checkout the consistency between the Driver and 
Navigation routines and to evaluate the larger quantization effects on 
the Navigation routine. Procedures for making these runs were described 
in Section 4.3.1, Reference 2. Expanded output word (0PW) lists for 
both runs and output for the first three seconds for the nominal run and 
the first six seconds for the alignment run of the present program are 
included in Appendix B. 

4.2.3 Checkout of Calibration Program 

Simulation runs shall be made as described in Paragraph 4.3.2, 
Reference 2. 

4.2.4 Alignment Filter Simulation Runs 

Prior to the alignment test phase (Section 6), simulated alignment 
runs, using the driver routine, shall be made to determine the alignment 
filter convergence time with the increased instrument quantization of 
the BB DDH. 

Figures 4-1 and 4-2 show the results of previous simulation runs 
for various accelerometer quantization values, all smaller than the 
BB DDH quantization. 
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Figure 4-1. Azimuth Misalignment versus Time With No I MU Sway Motion 
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Figure 4-2. Azimuth Misalignment versus Time With High Amplitude, 
Low Frequency IMU Sway Motion 
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5. CALIBRATION TEST PHASE 


5.1 TEST OBJECTIVES 

The specific objectives of this test phase are to: 

• Validate the modified RSP calibration routine, the calibration 
procedures, and all sign conventions, with the BB DDH. 

• Obtain compensation coefficients for the RSP Align, Navigation, 
and Failure Detection Routines. 

• Obtain performance data on: 

BB DDH parameter stability off-nominal performance; 
hysteresis, position sensitivity, etc. 

Calibration method accuracy 

5.2 TEST REQUIREMENTS 

The specific requirements of this test phase, to satisfy the test 
objectives, are: 

• Validate compatibility of the calibration program, the BB DDH, 
and the test equipment interfaces. 

• Validate the sign conventions, term by term, in the measurement 
and computation of all compensation parameters. 

• Establish the accuracy of the calibration systems by 1) error 
analysis and 2) empirical confirmation. 

• Determine compensation parameter coefficients for main program. 
Verify consistency with prior values for the same instruments. 

• Obtain repeated and redundant data and analyze for: 

parameter repeatability 
measurement accuracy 
OA axis g sensitive drift 
position sensitivity 
hysteresis 

channel intercoupling 
temperature sensitivities 
other anomalies 
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• Determine the required recalibration interval for the balance of 
the test program. 

• Develop operator skill in performing calibrations (a significant 
factor in inertial equipment testing). 

5.3 TEST CONSTRAINTS 

5.3.1 Software Status 

The modified calibration program shall have been verified with test 
cases as defined in Section 4. 

5.3.2 Hardware Status 

The hardware shall have been verified to conform to the configura- 
tion of Section 2 and Reference 2, Section 3. 

5.3.3 Inputs 

The calibration test stimuli shall consist of earths rotation and 
gravitational acceleration in a six-position static test, and these 
plus rate table rotation in a three-position rate test. 

5.4 TEST DESCRIPTION 

Note: General procedures and test details will be contained in the 

Detailed Test Plan. 

5.4.1 Calibration Test Setup 

The nine required test positions, shown in Figure 5-1, are obtained 
by use of the adjustable test fixture described in Reference 2, 

Paragraph 2.3.3, and the two degrees of freedom of the GOERZ table 
(rotary axes and trunnion axis). Prior verification of the alignment 
(and stability) of the GOERZ rotary axis to the vertical in position 1, 
and to the horizontal and North in positions 4 through 9, is assumed. 
Trim of the BB DDH axes to the vertical and North after each remounting 
is accomplished with leveling adjustments on the BB DDH base and on the 
mounting fixture, while observing the alignment cube optically. Coarse 
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STATIC TEST POSITIONS 
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The RSP requires sense switch settings to indicate which position 
the DDH is in and six angles 0 (i,j,k) indicating the azimuth of the 
horizontal axes in the static test positions (4 through 9) Table 5-1 
shows the set of azimuth angles 0 (i,j,k). 

Table 5- 1 1 shows the correlation between sense switch settings, 
the CTFG (cal test flag), and test table positions in the STATIC TEST 
mode. Note that the measured sensor output subscripts 1 through 6 
(Figure 3-9, Volume II, Reference 2) agree with the CTFG bit numbers. 

Table 5-1. Azimuth of Assembly Axes 


Equation Symbol 

Nominal Azimuth Angle 

0 YXU 

90 degrees (EAST) 

0 YXD 

270 

0 ZYU 

0 

0 ZYD 

0 

G XAU 

90 

0 XZD 

270 


Table 5-II. Static Test Cal Flag (CTFG) Status vs Calibration Position 


CTFG Bit & Sense SW (1) 

Orientation 

Position Number Figure 5-^ 

1 

X up 

4 

2 

X down 

5 

3 

Y up 

6 

4 

Y down 

7 

5 

Z up 

8 

6 

Z down 

9 
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The gimbal of each SAR must be positioned prior to the static tests 
to locate the spin axis or output axes as shown in Table 5-1 II , and v 
Figure 5-2, for compatibility with the calibration program, 

^ I A A. -=* % vw** 

SAR Gimbal Alignment in Calibration 


Table 5-III. 


SAR 

DDH Reference Axis 

XI, X2 
Yl, Y2 
Zl, Z2 

Output axis along +Z 
Spin axis along +X 
Spin axis along -Y 



5.4.2 Calibration Sequence 
5. 4. 2.1 Nominal Sequence 

Sensor alignment angles to the BB DDH reference axes are measured 
in the rate tests, positions 1, 2, and 3. These values are required in 
computation of static test results. The presently programmed sequence 
is rate tests, then static tests. The nominal sequence is: 

Position 


1 . 

Rate Test Position 

1 

2. 

Remount Fixture 


3. 

Rate Test 

2 

4. 

Remount DDH 


5. 

Rate Test 

3 

6. 

Compute Rate Test Results 


7. 

Static Test 

4 

8. 

Static Test 

5 

9. 

Static Test 

6 

10. 

Static Test 

7 

11 . 

Remount DDH 


12. 

Static Test 

8 

13. 

Static Test 

9 

14. 

Compute Static Test Results 
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AXIS 


Figure 5-2 


BB DDH Sensor Orientation 





5. 4. 2. 2 Abbreviated Sequence 

A revised calibration sequence minimizing the number of remountings 
and the total test time, follows: 


Position 


1. Static Test 4 

2. Static Test 5 

3. Static Test 6 

4. Static Test 7 

5. Rate Test 3 

6. Remount DDH 

7. Rate Test 2 

8. Static Test 8 

9. Static Test 9 

10. Remount Fixture 

11 . Rate Test 1 


12. Compute Rate Test Results 

13. Compute Static Test Results 

This sequence requires considerable manipulation to exit from the 
static and rate test routines repeatedly; the calibration program logic 
could be changed to make this out of sequence operation much easier. 


The static tests, and computation, can be performed separately by 
inputting (from paper tape) previously measured gyro scale factors and 
direction cosines. 

5.4.3 Calibration Alignment Accuracy Requirements 


An alignment error analysis has not been performed (see 
Paragraph 5.5.1), but preliminary estimates of the desired calibration 
test alignment errors are: 


Alignment to North 
Alignment to Table Rotation Axis 
Table Rotational Axis to Vertical 
Error in Axis Inversion 


3 arc minutes 
15 arc seconds 
10 arc seconds 
10 arc seconds 
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5.4.4 Calibration Test Output Data 


The rate test and static test outputs are detailed in Tables 3-XV 
and 3-XVIII of Volume 1, Reference 3. They include: 


Rate Test 

Gyro Scale Factors (6) 

Gyro Direction Cosines (6x3) 

Gyro C Matrix (15x6) 

Static Test 

Accelerometer Scale Factor (6) 

Accelerometer Biases-4 each (24) 

Accelerometer Direction Cosines (6x3) 

Accelerometer C Matrix (bypassed) 

Gyro Bias (6) 

Gyro Input Axis Mass Unbalance (6) 

Gyro Spin Axis Mass Unbalance (6) 


The four computed bias values for each accelerometer are for X, Y, 
and Z vertical, and the average of all three. 

In addition, data is available for hand calculation of M gp (drift 
proportional to acceleration along the gyro output axis), and (if data 
is taken in positions 10 and 11 (intermediate to 8 and 9)), redundant 
determinations of some of the above parameters are available. 

5.4.5 Calibration Test Duration and Rates 

The rate tests are to be performed at 10 degrees/second, clockwise 
and counterclockwise. The static test duration shall be 240 seconds. 

The static test interval could be increased to reduce the gyro drift 
measurement uncertainty (.02 degrees/hour) due to sensor pulse quantiza- 
tion (4.9 seconds), but this accuracy is sufficient for the breadboard 
testing. 

A single calibration sequence, allowing time for shutdown, remounting, 
and temperature stabilization is estimated to require six hours. 
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5.4.6 Calibration Test Series 


It is anticipated that two to three weeks will be required to obtain 
sufficient data to satisfy all of the test requirements listed in 
Subsection 5.2. Testing will consist of repetitions of all or part of 
the calibration sequence described above, and special investigations as 
required. 

In the initial stages, raw data (pulse counts) will be read out 
frequently to confirm the calculations and particularly, to verify that 
the correct polarity is assigned to bias and g sensitive drift terms. 

The test duration should be kept flexible and should last until the 
requirements are satisfied. Probably no less than 10 successful 
calibrations will be performed. 

A flow diagram of the test and analysis sequence is shown in 
Figure 5-3. 

5.5 CALIBRATION TEST ANALYSIS AND EVALUATION 

5.5.1 Test Accuracy (Error) Analysis 

A measurement error analysis, combining statistically the 
uncertainties in timing, position, counting, quantization noise, etc., 
should be performed to establish the expected accuracy in deriving the 
sensor compensation parameters. 

5.5.2 Sign Validation 

Sign errors in scale factors and direction cosines will be obvious. 
Each bias term (and particularly the g sensitive terms) should be 
examined for polarity by (1) comparing results with prior instrument 
test data, and (2) examination of • the accumulated count data. Only if 
the error terms are very small will sign errors be undetectable. 

5.5.3 Parameter Stability Evaluation Stability 

Records of all calibration data should be maintained for analysis 
of repeatability and trending. 
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Figure 5-3. Calibration Test Flow Chart 



















Settling Time 


The warmup time required for stable operation after each remounting 
should be determined empirically. 

Anomalous Behavior 

Redundant measurements exist for accelerometer bias, and can be 
derived for gyro bias. This data should be examined for consistency 
and evidence of position sensitivity or Interchannel coupling. 

If settling time, hysteresis or thermal lag type errors appear or 
are suspected in the calibration data, then special measurements should 
be implemented, varying the order or direction of rotation, thermal 
profile, or sampling interval. The three week estimate for the calibra- 
tion test does not presume problems of this type. 
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6. ALIGNMENT TEST PHASE 


6.1 TEST OBJECTIVES 

In this test phase the accuracy and settling time of the self- 
alignment scheme is to be checked under static and (single axis) dynamic 
conditions. 

Further alignment evaluation, with instrument failures is contained 
in Section 8, FDDC Tests. Alignment evaluation in the van environment 
is described in Section 9, Van Navigation Tests. 

6.2 TEST REQUIREMENTS 

The specific test requirements are: 

• Validate compatibility of the alignment program and the BB DDH 
hardware. 

• Evaluate alignment settling time with the large BB DDH sensor 
quantization values. 

• Evaluate accuracy and settling time as a function of initial 
azimuth estimate error. 

• Evaluate alignment accuracy and settling time with single axis 
oscillatory motion representative of stationary aircraft sway. 

6.3 TEST CONSTRAINTS 

6.3.1 Alignment Software Checkout 

The modified driver program shall have been used to validate the 
alignment routine, and specifically to determine filter convergence 
(settling time) with the larger values of gyro and accelerometer quanti- 
zation. These simulation results will dictate the nominal alignment time 
to be used in laboratory and van testing. 

6.3.2 Valid Compensations 

Alignment tests shall be preceded by a (recent) calibration. Update 
of gyro bias is particularly critical to alignment accuracy. If the 
calibration test analyses show that normal calibration procedures and 
delays will result in a total bias uncertainty during alignment in excess 
of .05 degree/hour, special procedures should be used to reduce the 
uncertainty. 
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6.3.3 SAR Gimbal (Head) Position 


The standard SAR gimbal positions, to be established prior to each 
calibration run, are identified in Table 5-1 I I . 

6.4 TEST DESCRIPTION 

6.4.1 Static Tests 

6. 4. 1.1 Test Positions 

The alignment tests shall be performed on the GOERZ Table. The 
standard test position shall be X vertical (up), Y North, and Z West 
(position 1 in Figure 5-1). 

6. 4. 1.2 Test Duration 

The nominal alignment interval will be established from the results 
of simulation (Driver) alignment runs. 

6. 4. 1.3 Input Data 

The site dependent data and actual BB DDH (and table) orientation 
angles shall be input, as indicated in Reference 3, Volume 1, Table 3-XIX, 
and as modified in Reference 2. Runs will be made with the initial 
azimuth estimate, A^, in error by 0, 5, and 15 degrees. Compensation 
parameters shall be updated by the performance of a calibration. 

6. 4. 1.4 Output Data 

The principal real time output data will be the A matrix (body 
axis to VEN coordinates) which is updated by the alignment routine. 

6. 4. 1.5 Reruns 

Reruns, to be specified in the Detailed Test Plan, will be performed 
to evaluate: 

• Repeatability 

t Settling time versus initial azimuth error 

• Accuracy versus bias compensation error (a polarity check) 

§ Accuracy versus alignment attitude 


6-2 




6.4.2 Dynamic Tests 

With the GOERZ table trunnion at 90 degrees as in position 4, 

Figure 5-1, a low frequency horizontal axis oscillation, shall be used to 
simulate vehicle sway motion during alignment. The motion shall be 
terminated and the BB DDH returned to its starting position prior to the 
completion of alignment. The simulated sway motion shall be as specified 
in Table 6-1. BB DDH position, input and output data are as discussed 
under Static Tests. 

Table 6-1. Simulated Sway Motion 


Run 

Sway -Amplitude 
0-Peak 

Frequency 

H z 

1 

1 .6 degrees 

0.1 

2 

0.6 degrees 

0.25 

3 

0.3 degrees 

0.5 


6.4.3 Test Sequence 

Figure 6-1 is a flow diagram of alignment tests (without failures). 
It is anticipated that two weeks will be adequate time to evaluate the 
alignment routine and debug the procedures. 

6.5 ALIGNMENT TEST ANALYSIS AND EVALUATION 

The static and sway alignment data will be evaluated for consistency 
with simulation results under the various test conditions. 

Alignment errors (in VEN coordinates) will be computed from the 
alignment attitude reference matrix A (body axes to VEN) and the known 
body to VEN attitude matrix K, by the cross product of the unit vectors: 
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Accuracy 
Settling Time 


* For various values of initial 
azimuth error. 


Figure 6-1. Alignment Test Flow Chart 
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Similarly, alignment errors about body axes 
correlation with sensor errors. 



can be computed, for 


The first of these computations is implemented in the existing pro- 
gram but a minor modification is required to make it usable in the 
real time align mode. 


Convergence time can be monitored on the real time printout by 
noting when alignment is essentially complete, e.g., A matrix indicates 
alignment within .5 degrees of final value. 


Alignment program output data will be recorded on magnetic tape. 
A useful external data evaluation tool would be a computer plotting 
routine for alignment errors. 
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7. NAVIGATION TEST PHASE 


7.1 TEST OBJECTIVES 

The objectives of this test phase are to validate the compatibility 
of the navigation program and the BB DDH hardware and to evaluate 
navigation performance without instrument failures. 

7.2 TEST REQUIREMENTS 

The specific test requirements are: 

§ Validate hardware/software compatibility and compensation 
adequacy. 

• Evaluate navigation accuracy for static and single axis rota- 
tional tests, compare with expected values. 

• Establish baseline data for FDDC tests. 

7.3 TEST CONSTRAINTS 

7.3.1 Sensor Compensation 

The calibration interval required to maintain valid compensation 
coefficients during this test phase will be determined from analysis 
of the calibration test phase data analysis. 

Navigation accuracy is particularly dependent upon accurate com- 
pensation for horizontal gyro drift. Total bias uncertainty (compen- 
sation error) in the navigation -test orientation, should not exceed 
.05 degrees/hour per axis. Special calibration procedures should be 
used if stability between normal calibrations is inadequate. 

7.3.2 A1 ignment 

Navigation tests will normally be preceded by self-alignment. In 
special cases, designated in the Detail Test Plan, the alignment mode 
may be bypassed. The navigation program will then initialize to the 
input values which must accurately represent the BB DDH orientation. 

7.3.3 Gimbal Positions 

The standard SAR gimbal positions, to be established prior to each 
navigation run, are identified in Table 5-1 1 1 . 
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7.4 TEST DESCRIPTION 


Note: Effective test analysis and error identification .requires 

that very simple navigation profiles be employed. 

7.4.1 Test Position 

The standard static test position shall be X up, Y north, and Z 
west (position 1 in Figure 5-1). 

7.4.2 Test Duration 

Eighty-four minutes are required for one cycle of errors related to 
gravity. Shofter term tests are useful for examination of certain 
effects. The two principal navigation intervals have, therefore, been 
chosen as 20 minutes and 90 minutes. 

7.4.3 Input Data 

The BB DDH (test table) orientation must be input. If different 
from the nominal program values. Nominal program values must have been 
established for the site geodetic parameters and other constants listed 
in Table 3-XXII, Volume I, Reference 3. Valid compensation parameters 
must be in the program from a recent calibration. 

7.4.4 Test Method 

For each navigation run (prior alignment assumed): 

• insert initial attitude 
t insert run duration 

• initiate navigation 

t introduce table motion (if required) 
t observe velocity, position, or attitude errors 

• tape selected word list 

• terminate table motion 

• terminate navigation 

• printout final errors 
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7.4.5 Output Data 

Real time output data is listed in Table 7-1. 

Table 7-1. Navigation Test Real Time Output Data 


Sense 
Swi tches 

Word 1 

Word 2 

Word 3 

Parameter, 

VEN Coordinates 

3 

4 





ON 

OFF 

Vert 

East 

North 

Velocity Errors 

OFF 

ON 

East 

North 

Altitude 

Position Errors 

ON 

ON 

0 V 

0 E 

0 N 



Taped output data selected from the major cycle output word list, 
Table 3-XXIII, Volume 1, Reference 3, will include: 

• Latitude, longitude, VEN position, and velocity 

• VEN position velocity and attitude errors 

• Attitude matrix 

• Test table rotation pulse accumulation 
Minor cycle taped data will include: 

t 6 gyro pulse accumulations 

• 3 accelerometer pulse accumulations 
t table angle 

• discrete word 

Taped data may be used in the tape input mode to make repeated 
navigation runs with the same sensor data. 

7.4.6 Navigation Test Accuracy 

Navigation accuracy will be limited by the instrument performance In 
the BB DDH thermal and mechanical environment. There is no closed loop 
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temperature control. Lack of compensation for the SAR gimbal rotation 
will add some error but this will probably be masked by normal SAR 
Instability. 

Based on estimates in the error analysis section, the expected 
system navigation performance is in the 5 to 15 kilometers/hour range, 
depending on the bias update accuracy. These estimates are for simple 
navigation profiles, do not include dynamic errors, and presume constant 
ambient temperature. 

7.4.7 Test Sequence 

The test variables (to be specified in the Detailed Test Plan) 
will include: 


Parameter 


Nominal Value 


Test Duration 
Orientation(s) 

Rates 

Alignment Method (Self or External) 
Altitude Damping Coefficients 


90 minutes 
Position 1 
10 degrees/second 
Self 



Less time should be required for this test phase than for the align- 
ment test phase as the main program/BB DDH interfaces will have been 
debugged. 

7.4. 7.1 Static Tests 

Fixed attitude tests will be performed first to verify proper 
operation and compensation. Enough runs will be made to obtain statistics 
on repeatability. 

7. 4. 7. 2 Dynamic Tests 

The basic dynamic test will involve a single rotation. This will 
demonstrate the errors to be encountered in van testing, but is still 
simple enough that reasonably accurate error predictions can be made. 
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The test table motion and rates will be specified in the 
Detailed Test Plan. A typical test profile will be: 

1) Initiate navigation 

2) Rotate 90 degrees at 10 degrees/second 

3) Navigate in this fixed position for 90 minutes, terminate 

A series of tests for different rotation angles will be made. 

The second category of dynamic tests shall be static navigation with 
angular vibration about a horizontal axis, to evaluate system operation 
in the presence of simulated road motion. The vibratory amplitude shall 
be 1 degree/second at 5 Hz, or at the maximum frequency of which the 
GOERZ table is capable. When (If) actual measurement data on the van 
angular rotation environment becomes available, the laboratory test rates 
should be revised accordingly. 

Figure 7-1 is a navigation test flow chart. 

7.5 TEST ANALYSIS AND EVALUATION 

Time profiles of the VEN navigation errors listed in Table 7-1 will 
be the primary tool for navigation analysis. Computer plotting routines 
for these parameters (from tape) should be developed. 

For each run, accuracy predictions should be made as in the error 
analysis. Appendix A, Section A4.‘ By use of the error propagation 
equations and plots of Appendix A, identification of major sources of 
error in navigation runs should be possible. 

Tapes from the runs make it possible to repeat the input data if a 
rerun is desired. This will be particularly useful in evaluating van data 
and FDDC performance. 
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8. FDDC TEST PHASE 


8.1 TEST OBJECTIVES 

The objective of this test phase is to evaluate performance of the 
FDDC logic with actual sensor outputs. 

8.2 TEST REQUIREMENTS 

The specific requirements of this test phase are: 

• Establish empirically the minimum filter constants required to 
prevent false alarms. 

• Evaluate failure detection time versus predicted time in the 
presence of noise. 

• Evaluate navigation error with sensor failures. 

• Evaluate navigation error with unfailed, but degraded sensors 
(below the detection threshold). 

• Evaluate multiple instrument failures. 

t Demonstrate internal monitor failure logic capability in isolating 
a third failure. 

• Evaluate FDDC logic susceptibility to instrument shock. 

8.3 TEST CONSTRAINTS 

8.3.1 Sensor Compensations 

Valid sensor compensation constants shall be in the program prior to 
each run. See Paragraph 7.3.1 

8.3.2 Filter Constants 

Initial values of the gyro filter constants (Kfl g and Kf2 g ) have been 
established, based on the model of expected sensor performance and noise. 
See Paragraph 8.4.3. 

8.4 TEST DESCRIPTION 
8.4.1 General 

In general, the FDDC tests will consist of repeats of the prior align 
ment and calibration tests (static and dynamic) with the introduction of 
one or more sensor failures. The effects of catastrophic failures are 
evaluated first, then the effects of degraded sensor performance. 
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8.4.2 Test Duration 


Duration will be specified in the Detailed Test Plan. Some tests 
will be terminated as soon as a failure is detected. 

8.4.3 Input Data 

As for normal alignment and navigation tests, with the addition of 

the prefilter constants Kfl (time constant) and Kf2 (gain or threshold). 

y y 

Normal values (from Appendix C) are: 

Kfl = .99992 
9 

Kf2g = 80 seconds/rad 

8.4.4 Test Method 

For each real time test run: 

• Set filter constants. 

• Initiate alignment or calibration run as in Section 6 or 7. 

• Introduce table motion (in dynamic tests). 

• Tape record FDDC and navigation data. 

• Introduce sensor failures. 

• Observe the gyro failure state word (FSG) for time of failure. 

• terminate motion, terminate run. 

For taped data reruns: 

• Set filter constants. 

• Start tape. 

• Initiate align or navigation run. 

• Terminate run, terminate tape. 

8.4.5 Output Data 

The gyro failure state word can be monitored real time. The special 
output data which will be taped is listed in Table 3-XXVII, Volume 1, 
Reference 3. 
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Minor cycle (sensor 40 millisecond accumulations) will be tape 
recorded, so that repetitions of the same run can be made for filter 
parameter evaluation. 

8.4.6 Alignment Test Sequence 

The alignment tests will be repeated with catastrophic sensor errors, 
using output monitoring only. The sequence is: 

• Static 

Repeat static alignment static test without failure as per 
Paragraph 6.4.1 (-Z axis east). Repeat static test with Z1 gyro 
failed at 30 seconds. (Disconnect gyro output from IFE). 

Repeat with Z1 failed at 30 seconds, Z2 failed at 35 seconds. 

• Dynami c 

Repeat dynamic alignment test, run 2 of Paragraph 6.4.2, with 
no sensor failure. Repeat dynamic alignment one (easterly) gyro 
failed at 30 seconds. Repeat with one gyro failed at 
30 seconds, a second failed at 35 seconds (both gyros will be 
those nearest the east pointing DDH axis). Test evaluation is 
discussed in Paragraph 8.4.8 below. 

8.4.7 Navigation Test Sequence 
The general sequence is: 

t Verification of nominal filter threshold (gain) and time 
constants (no false alarms). 

• Catastrophic failures (1 and 2 gyros failed) output comparison 
only: 

• Static navigation test 

• Dynamic navigation test (subparagraph 7. 4. 7. 2). 

• Catastrophic failures (1; 2, and 3 gyros failed), using internal 
monitoring only, static test only. 

• Catastrophic failures (3 gyros failed), using both output and 
internal monitoring. 

• Static test only 

• Soft failures, using output comparison only, varying filter 
parameters and sensor error level 

• Shock susceptibility 
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The time of failure, the channels to be failed, and the magnitude 
(drift rate) of soft failures will be specified in the Detailed Test Plan 
The BB DDH control panel is capable of introducing soft failures of 
.4 and 4 degrees per hour in any channel. 

Figure 8-1 is a flow diagram of FDDC Testing. 

8. 4. 7.1 Threshold Verification 

The first test sequence will establish the margin of the nominal 
failure threshold level by lowering it until false alarms occur due to 
normal system performance. The nominal threshold will be raised if 
necessary to eliminate false alarms during static operation and during 
a turn of 90 degrees at 10 degrees/second. 

8. 4. 7. 2 Soft Failure Testing 

This is the test phase where parametric experimentation will occur. 
The time required to detect a failure will be measured as a function of 
(1) instrument error and (2) failure threshold setting. The magnitude of 
navigation errors accumulated before detection of a failure will be 
measured. Probable values of instrument error and threshold setting are: 

Instrument Error 0, .4, and 4.0 degrees/hour 

Threshold Setting .5, 1, and 2 times nominal 

To provide repeatable data for evaluation of filter parameter 
changes, the input data shall be from tape, except for the first run. 

8. 4. 7. 3 Shock Susceptibility 

The RSP failure detection logic locks up any of the 15 elements of 
the V vector (test signals) on the occurrence of a single check in 
excess of the threshold. There is some probability, therefore, of an 
accumulation of non-zero elements due to transients or shocks to the 
sensors, and ultimately switching out a sensor which had not permanently 
failed. A test shall be performed to examine the failure states 
(15 element V vector) with repeated shocks to the DDH, to evaluate the 
vulnerability of this FDDC logic scheme to noise and transients. 
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Figure 8-1. FDDC Test Flow Chart 










8.5 TEST ANALYSIS AND EVALUATION 

8.5.1 Selection of Nominal Filter Constants 

Nominal initial values for the FDDC prefilter gain and time constant, 
derived in Appendix C, are: 

Kfl = .99992 
9 

Kf2 = 80 
9 

The corresponding time constant is 500 seconds. The detection 
threshold on the test signals, V. , is .2 degrees/hour. Revisions to the 
nominal values may be based on the observed gyro (compensated) errors at 
the time of the FDDC test, using the method of Appendix C. 

8.5.2 N avigation Accuracy With Sensor Failures 

Predictions shall be made for navigation errors with 1, 2, or 3 
sensors failed, using the methods of Appendix A, and results of tests 
shall be compared with predictions as discussed in Subsection 7.5. 
Anomalies shall be investigated by identifying the errors with the 
characteristic error functions given in Appendix A. 

A special test and analysis that shall be performed to verify 
analytical results, is to increase sensor (gyro bias) errors to the 
maximum values that will not trigger the failure detection logic, then 
measure the navigation system errors in this degraded mode. The 
Failure Simulator routine will be used to inject these errors. 

Failure detection times versus sensor errors and threshold levels 
should be evaluated against expected values using the methods of 
Appendix C. 

From Appendix C it can be seen that for times long compared to the 
filter time constants (>30 minutes) the navigation error is limited 
(approximately) to the detection threshold rates, or .2 degrees/hour 
(approximately 22 Km/hour). 
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For times short relative to the time constant (<3 minutes) the 
navigation error is limited to the angle e = x/Kf = .028 degrees or 
approximately 3.2 Km 


9. VAN NAVIGATION TEST PHASE 


9.1 TEST OBJECTIVES 

The objective of this test phase is to demonstrate the system align- 
ment, navigation, and failure detection/correction capability in a 
moving vehicle. 

9.2 TEST REQUIREMENTS 

The specific test requirements 

t Establish an operating land navigation system with alignment and 
initialization techniques, in a mobile van. 

• Obtain taped sensor data from typical vehicle test runs for 
repeated use in navigation program reruns in the laboratory 
computer. 

• Evaluate self -alignment capability in the mobile van. 

• Evaluate navigation capability in the van. 

• Evaluate FDDC operation in the mobile van environment. 

9.3 TEST CONSTRAINTS 

9.3.1 Configuration 

Real time navigation will not be done since the navigation computer 
will not be carried in the vehicle. Taped data will be used for input 
to the DDP-124. 

Thermal and dynamic environment will not be tightly controlled. No 
privision now exists for measurement of angular vibration spectra in the 
van. 

9.3.2 Sensor Compensation 

Accuracy of the gyro compensation will probably be the limiting 
factor in navigation accuracy. 

Prior to navigation runs if the BB DDH attitude is stationary and 
precisely known, by external alignment, total sensor biases can be 
extracted from the taped data, and used for compensation updates prior to 
making the taped navigation run. Procedures for determining gyro bias in 
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this manner have not been established. The accuracy will be determined 
by the knowledge of BB DDH attitude at the time the reference data is 
recorded. 

9.3.3 Gimbal Positions 

Prior to vehicle alignment or navigation runs, the SAR gimbals will 
be set to standard reference positions (to be specified in the 
Detailed Test Plan). 

On-board counters will be used to monitor and set the gimbal 
positions . 

9.4 TEST CONFIGURATION 

9.4.1 BB DDH Installation 

The BB DDH will be mounted in the van with the X axis nominally 
vertical and the Y or Z axis parallel to the longitudinal (fore-aft) 
axis of the vehicle. Vibration isolators will not be used. 

Provision will be made for optical access to the BB DDH alignment 
prism for external azimuth determination. 

9.4.2 Van Instrumentation 

The van instrumentation for measurement of auxilliary test data 
will consist of: 

t tape recorder (14 track) 

t gimbal angle readout (6) 

• speedometer 

• odometer 

• thermometer 

• vibration monitors 

Considering the expected system accuracy (5 to 15 Km/hour) the 
vehicle odometer and speedometer will provide adequate reference data. 

A fifth wheel type precision odometer/speedometer is not required. 
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Accelerometers for linear vibration measurement are provided on the 
BB DDH. No provision currently exists for angular vibration measurement, 
although this data is potentially more useful (in the event of 
significant coning or angular rectification type errors). 

9.5 TEST DESCRIPTION 

9.5.1 Test Duration 

The basic test interval shall be 90 minutes, to encompass one 
Schuler period. By establishing frequent checkpoints, the recorded data 
can also be used for evaluating navigation for shorter intervals. The 
necessary alignment interval will be determined from simulations 
(Paragraph 4.2.4). 

9.5.2 Navigation Profiles 

Error analysis will be enhanced by simple vehicle routines, e.g., 
right angle control led-rate turns, and flat terrain. Some runs will be 
made with no maneuver constraints (other than maximum turning rates) but 
only gross evaluations will be possible. The course will be arranged to 
provide frequent surveyed checkpoint crossings. 

Vehicle velocity will only be limited by the dynamic (vibration and 
shock) environment. Thirty mph (forty-eight Km/hour) will be the nominal 
cruising speed. 

Vehicle turning rates will normally be limited to 90 degrees at 
10 degrees/second to conform to laboratory tests and to avoid triggering 
the failure detection logic. In any case the turning rate must be less 
than 1 radian/second (the SAR limit) to avoid loss of inertial reference. 

9.5.3 Bias Update 

If an on-board bias update is desired prior to self-alignment or 
navigation, sensor data should be recorded for approximately 10 minutes 
in a stable precisely known attitude. A combination of BB DDH bubble 
levels, optical azimuth measurement, and vehicle maneuvering will be 
required to establish the attitude. 
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9.5.4 Navigation Coordinate System Alignment 

Self-alignment will normally be used. However, by inputting the 
known BB DDH attitude as determined in Paragraph 9.5.3, if the self- 
alignment routine is bypassed, the navigation program will initialize 
(align) to the externally measured coordinates. 

The externally measured alignment data will be used as a reference 
in evaluating self-alignment accuracy. 

9.5.5 Input Data 

When making computer navigation runs with taped data, the site 
dependent data shall be loaded as for laboratory navigation tests, 
Paragraph 7.4.3. 

The input data normally representing BB DDH orientation on the GOERZ 
E T L 

(£, Rj^, R l , and Rj) can be manipulated to represent the inertial 
orientation of the BB DDH, as measured optically (Paragraph 9.5.3) at 
the beginning of the navigation run. These inputs can be bypassed, as 
can external alignment and bias update, in van testing, with some decrease 
in accuracy. 

9.5.6 Output Data 
9. 5. 6.1 Van Data 

Tape recorded data will include: 
t 6 SAR outputs 

• 3 accelerometer outputs 

• time of day 
t clock 

• vibration instrumentation 

Continuous tapes are desired. If necessary, tape changes can be made 
while stationary. 

Manually recorded data will include: 

• event-time log 

• checkpoint ID, time 
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• velocity 

• odometer heading 

• gimbal angle readout 

• temperature 

9. 5. 6. 2 Navigation Program Output Data 

This will be the same as for laboratory tests. Paragraph 7.4.5. 
9.5.7 Test Method 

9. 5. 7.1 Van Navigation Recording Run 
For each run: 

• allow system temperature stabilization 

• initialize SAR gimbal angles 

• align vehicle with bubbles, external optics 

• record sensor data for alignment interval (to be determined) 
avoiding all motion 

• realign gimbals 

• drive prescribed course, recording sensor data 

• induce sensor failures or errors at prescribed time 

• stop at checkpoints and for tape changes 

• record environmental and event data 

• • * 

• return to starting point and realign 

• record static sensor data 

9 . 5 . 7 . 2 Laboratory Computer Runs 

• Initialize alignment program. 

• Insert valid compensation coefficients. 

• Run data tape (from van). 

* 

Optional 

Failures (hard and soft) will generally be introduced at playback. 
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• From gyro sensor accumulations and known BB DDH alignment, compute 
gyro bias compensation error*. 

• Update gyro bias compensation. 

• Enter alignment routine. 

• Record alignment data. 

• Enter navigation routine. 

• Induce failures as specified. 

• Record data. 

• Terminate run. 

9.5.8 Test Sequence 

9. 5. 8.1 General 

After the van operation, procedures, recording equipment are 
debugged, relatively few navigation runs need be recorded. These runs 
(or portions of them) can be played back repeatedly, introducing any type 
of failure at will, to evaluate the align, navigation and FDDC routines. 

In these runs, the compensation parameters, alignment method, and filter 
parameters can be varied. 

9. 5. 8. 2 Static Runs 

The first recorded van data will be in a stationary position with 
known BB DDH attitude for 90 minutes. This will be used to perform 
alignment and static navigation runs for comparison with laboratory test 
results. 

9. 5. 8. 3 Comparison Run for Laboratory Dynamic Tests 

A run with a single 90-day rotation, preceded and followed by a 
static (alignment interval) should be performed for comparison with 
laboratory dynamic tests. 

9. 5. 8. 4 Maneuver Series 

Tapes shall be recorded for a series of runs progressing in 
complexity from straight line constant speed to multiple turn variable 
★ 

Optional 
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speed. On some of these runs, internal monitor and gyro bias type 
failures should be introduced. Each run should be preceded and followed 
by a static data collection interval (for alignment or, if optically 
aligned, bias determination). 

Ten or fewer 90-minute navigation tapes should provide all of the 
required data for navigation evaluation, except for the investigation 
of anomalies. 

9. 5. 8. 5 Laboratory Computer Runs (Alignment and Navigation) 

The evaluation sequence shall be: 

• static alignment 

• static alignment with failures 

• static navigation 

• straight line navigation 

• maneuvering navigation 

• maneuvering navigation with failures 

• FDDC tests 

The FDDC tests will establish usable filter constants, and the 
failure detection effectiveness under field conditions. 

Figure 9-1 is a flow diagram of the van test phase showing van 
tests and laboratory navigation runs. 

If tape changes are necessary in navigation runs, a special procedure 
for suspending the navigation program in the absence of inputs will be 
required. 

9.6 TEST ANALYSIS AND EVALUATION 
9.6.1 Methods 

The analysis methods will be the same as for laboratory alignment 
and navigation tests. 

Accuracy, settling times, and failure detection thresholds and 
effectiveness will be re-evaluated in an operational environment. 
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Figure 9-1. Van Test Phase Flow Chart 





















Computer plotting of navigation data will be especially valuable in 
attempting to relate maneuvers, navigation errors, and their sources. 

9.6.2 Accuracy 

The expected accuracy is less than in laboratory tests, probably in 
excess of 10 Km/hour for unconstrained maneuvers. 
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APPENDIX A 


ERROR ANALYSIS 

REDUNDANT SENSOR STRAPDOWN IMU TESTS 


A1.0 INTRODUCTION 

The redundant sensor strapdown IMU tests will be performed in a test 
laboratory with a DDP-124 computer and in a mobile van. The test article 
will be the MFSC Breadboard Dodecahedron (BB DDH). The tests will be 
performed with the Redundant Sensor Computer Program (RSP) in the DDP-124. 

The RSP was prepared for a series of tests which were to be per- 
formed at ERC with similar equipment. There are several differences 
between the equipment that will be used and the corresponding ERC equip- 
ment that will affect system performance. 

• The BB DDH contains single axis references (SARs). 

The ERC test article contained single axis 
strapdown gyros. The SAR gyro is free to rotate 
with respect to the DDH structure, and the input 
axis mass unbalance drift and the cross coupling 
caused by SAR internal misalignments are both 

a function of SAR gimbal angle. 

• The BB DDH contains three accelerometers with 
input axes nominally parallel to the DDH body axes. 

The ERC unit contained six accelerometers with 
input axes nominally parallel to the gyro input 
axes . 

• The quantization of the BB DDH accelerometers and 
SARs is much greater than that of the sensors in 
the ERC unit. The quantization of the BB DDH accel- 
erometers and SARs is 0.4 f ps/pul se and 4.9439 

arc seconds/pulse, respectively. (The quantization 
of the accelerometers will be reduced by a factor 
5 to 10). The quantization of the corresponding 
sensors in the ERC unit is 0.001286 fps/pulse and 
0.05 arc sec/pulse. 
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The laboratory tests will be performed with the GOERZ model 500 
Test Table that was to be used in the ERC tests. In the forthcoming 
tests, the two axis fixture that was planned for the ERC tests will not 
be used. The test article will be mounted with a fixture such as that 
described in Section 2. 3. 3. 2 of Reference 2. This type of fixture 
requires remounting during each calibration. 

There are several characteristics of the RSP that will affect 
performance. 

• The alignment filter in the RSP was designed for 
the quantization of the ERC unit. Convergence time 
will be increased by the larger quantization of the 
BB DDH. 

• Gyro drift compensation in the RSP was designed 

for the ERC unit. There is no provision for applying 
gyro input axis mass unbalance compensation and SAR 
internal misalignment compensation as a function 
of SAR gimbal angle. 

• There is no provision in the RSP for updating the 
gyro drift compensation during prenavigation 
alignment. 

t There is no provision in the RSP for filtering the 
gyro and accelerometer outputs during calibration. 

The performance capability of the inertial navigation system may be 
improved by suitable RSP changes. Such changes, however, are undesirable 
and will be held to a minimum. If changes can be avoided by test 
design, this approach is preferred. 

In the interest of minimizing RSP changes, a navigation accuracy 
goal of 5 n.mi./hr will be assumed for the current test series. Cali- 
bration and alignment accuracies will be consistent with this goal. If 
greater accuracy should be required in later tests, additional RSP 
changes may be required at that time. 



A brief analysis of the principal errors that will be generated in 
the current series of tests is presented in this appendix. The purposes 
of this analysis are as follows: 

• Identify the RSP changes that will be required 
because of performance considerations. 

• Provide test design guide lines. Define constraints 
that will avoid RSP changes. Provide a basis for 
determining calibration, alignment and mounting 
fixture accuracy requirements. 

• Provide data that may be used in the evaluation 
of test results. 

The analysis in this appendix is necessarily limited. This 
analysis should be extended by additional analyses as the program pro- 
gresses. 

A2.0 NAVIGATION ERROR ANALYSIS 
A2 . 1 Introduction 

Navigation errors will first be discussed for a nonredundant three 
axis strapdown IMU. It will be assumed that the IMU is oriented with 
two axes (Y and Z axes) level and the third axis (X axis) vertical. 

With the exception of the attitude drift that results from errors in the 
transformation of the earth rate vector from body coordinates to inertial 
coordinates, a nonrotating earth will be assumed. 

Once the navigation errors of a nonredundant IMU have been developed, 
the discussion will be extended to a dodecahedron with no failures. The 
performance degradation associated with fault detection, diagnosis and 
correction (FDDC) will be discussed in Section A3.0. 

A2.2 Gyro Drift, Level Gyros 

If the orientation of the IMU is fixed, uncompensated constant 
drift of one of the level gyros will cause the following attitude error: 

E, 

(A -’) 
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where: E, = attitude error, rad 

<p 

E-j = level gyro drift, rad/sec 
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R 
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2 

gravitational acceleration, ft/sec 
earth radius, ft 
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The attitude error E^ thus is periodic with period T Q (84 min). 

With the orientation of the IMU fixed, uncompensated drift of one 
of the level gyros will result in the following velocity and position 
errors: 

E y = E-|R[i - cos V ] (A-4) 


E = E, R ft - — sin u> t| 
p I L u)„ o J 
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where: E y = velocity error, ft/sec 
Ep = position error, ft 

Y gyro drift will result in velocity and position errors in the Z axis. 
Conversely, Z gyro drift will result in velocity and position errors in 
the Y axis. 
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The velocity and position errors both contain periodic terms with 
period T . The velocity error also contains the constant term E-jR. 

The corresponding position error term is the linear term E-|Rt. 

The error sensitivity of the constant term in Equation (A-4) ex- 
pressed in units of (n.mi ./hr)/(degree/hour) is as follows: 


(2.09x10 7 )(4.85x10~ 6 )(3,600) 

6,000 


= 60.8 


n.mi ./hr 
deg/hr 
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where: E 1 = velocity error, n.mi. /hr 

E-, 1 = level gyro drift, degree/hour 

The attitude, velocity and position errors (assuming no IMU orien- 
tation change) for a gyro drift of .01 degree/hour are plotted in 
Figures 3-1, 3-2 and 3-3 of the main body of this document. 

A2.3 Azimuth Alignment Error 

In the RSP, the direction cosine matrix relating inertial coordinates 
to body coordinates (A matrix) is initialized by means of the alignment 
subroutine. An error in the initialization process can be viewed as an 
initial alignment error. Such an error may cause an error in the trans- 
formation of the earth rate vector from body coordinates to inertial 
coordinates and thus cause attitude drift. 

The attitude drift in body coordinates caused by an initial azimuth 
alignment error is as follows: 

E^ = E^ ft cos A sin p (A-7) 


where: 


E^ = attitude drift, rad/sec 

E^ = azimuth alignment error, rad 

ft = earth's rotation rate, rad/sec 
A = test latitude, deg 

4 * = angular displacement of level body axis from north, deg. 



The earth rate transformation error will cause attitude, velocity 
and position errors. These may be computed with Equations A1 , A4 and A5. 

The error sensitivity of the constant velocity error term (Equa- 
tion A4) expressed in units of (n.mi ./hr)/arc seconds is as follows: 

(x is assumed to be 35 deg and ip is assumed to be 90 deg) 

E v ' = R fi cos 35° sin 90° 

r ' - (2.09xl0 7 ) (4.85xlQ" 6 ) ( 3,600) (1 5) ( 4.85x1 0~ 6 ) ( .819) 

t 2 “ 6,000 

= .0036 n - - mi ‘/ hr 
arc sec 

where: E 2 ' = azimuth alignment error, arc second. 

In the RSP, the reference coordinate system is aligned in azimuth 
by gyrocompass ing. Uncompensated drift of one or both of the level gyros 
will thus cause an azimuth alignment error. This alignment error, in turn, 
will cause an earth rate transformation error that will at least partially 
offset the uncompensated gyro drift. Once navigation is initiated, this 
error offset will continue until the orientation of the IMU is changed. 

A change in IMU orientation will cause the earth rate transformation 
error transformed to body coordinates to change, and thus cause the total 
drift to increase. 

This error mechanism will be illustrated by a special case. Assume 
that during alignment and during the initial portion of the navigation 
period, the IMU is positioned with Z axis north and Y axis east as shown 
in Figure A-l . At the termination, of alignment, the drift in each body 
axis will be: 


= E.| - E 2 n cos x sin ip 
where: E = total drift, rad/sec 
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In the Y axis (</i = 90 deg). 


E y = E ly “ E 2 n cos X 

The alignment process will result in a misalignment E 2 such that 


and 


In the Z axis (ip = 0 deg). 


E 2 fi cos x = E-|y 


E y = 0 


E z * E iz 


Assume that during the navigation period, the IMU is rotated 90 deg 
about the vertical axis as shown in Figure A-2. In the Y body axis after 
rotation = 180 deg). 


E 


y 



In the Z axis = 90 deg). 



The 90 deg rotation has thus caused the magnitude of both E y and E z to 
increase. (E^ and E-j z are assumed to be statistically independent). 

In addition to the foregoing errors, an azimuth alignment error 
will cause a position error that is a function of distance from the 
starting point. 

E p = E 2 d (A-9) 

where: d = distance from starting point, ft. 

This is the total error (the rss of the Y and Z axis errors). This error 
will be much smaller than the position error resulting from the earth 
rate transformation error. 
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A2.4 Gyro Drift, Vertical Gyro 

Uncompensated constant drift of the vertical gyro will cause an earth 
rate transformation error that changes linearly with time. The error 
component in each level body axis will be as follows: 


E r = E^t Q cos x sin . (A-10) 

where: E^ = vertical gyro drift, rad/sec. 

If the orientation of the IMU is not changed, the attitude, 
velocity and position errors that result from the foregoing error are as 
follows: 

E., q cos x sin ip r 

= — 5 1 - cos u> t 

<J> 2 I o 

w o 
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in t sin w t 

r U) 0 
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E = E^, R a cos x sin tj> 

P O 


| ^2 P ' cos “ 0 t) 

CO 
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The foregoing equations apply to both the Y and Z axes. For \p defined 
for the Y axis, E applies to the Y axis, and E v and E p apply to the 
Z axis. Conversely, for ip defined for the Z axis, E^ applies to the 
Z axis, and E v and E p apply to the Y axis. 

The attitude, velocity and position errors all contain periodic 
terms with period T Q . In addition, the attitude and position errors 
contain constant terms, the velocity error a linear term, and the position 
error a second order term. 

The error sensitivity of the linear term in Equation A- 1 2 evaluated 
at t = 45 min and expressed in units of (n.mi ./hr)/(degrees/hour)is as 
follows: (x is assumed to be 34 deg and ip is assumed to be 90 deg) 



= R n (cos 35°) (sin 90°)(45)(60) 



t = 45 min 

(2. 09x1 0 7 ) (4.85x1 0~ 6 ) (3 .600) (15) (4.85x1 0~ 6 ) ( .819) (45) (60) 

6,000 


q 7 n.mi ./hr 
' deg/nr 


(A-14) 


where: Eg 1 = vertical gyro drift, degree/hour 

The attitude, velocity and position errors (assuming no IMU orien- 
tation changes) for a vertical gyro drift of .01 degrees/hour are plotted 
in Figures 3-4, 3-5 and 3-6. 

In addition to the foregoing errors, vertical gyro drift will cause 
a position error that is a function of velocity history during the navi- 
gation period. For the special case of IMU velocity (magnitude and 
direction) constant, the position error is as follows: 


F H + 



(A-15) 


This is the total error (the rss of the Y and Z axis errors). This error 
will be much smaller than the position error resulting from the earth 
rate transformation error. 

A2.5 Gyro Misalignment 

If the IMU is rotated about the X axis, a misalignment of the Y gyro 
in the XY plane will cause an attitude error in the Y gyro axis. The 
same will be true for any other axis pair. Thus, 

S = E 4 e (A - 16) 

where: = gyro misalignment, rad 

e = IMU angular displacement, rad. 

Rotation about the vertical axis will cause a level error in each of 
the two level gyro axes. The attitude, velocity and position errors 
resulting from a level error are discussed in Section A2.7. 
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Rotation about one of the level gyro axes will cause both an azimuth 
error and a level error. The attitude, velocity and position errors 
resulting from an azimuth error are discussed in Section A2.3. 

The error sensitivity for rotation about the vertical axis expressed 
in units of n. mi. /arc seconds is as follows: (A rotation of 360 deg 

is assumed.) 

E 1 

- Re 

4 7 6 

(2.09x10] (4.85x10 d )(2tQ 

6,000 

= jq6 (A-17) 

arc sec 


where : 



= position error, n.mi. 

= level gyro misalignment, arc sec 


In addition to the foregoing errors, misalignment of the gyros with 
respect to the accelerometers and misalignment of the level gyros with 
respect to each other will cause the gyros to sense an erroneous earth 
rate component. This is a consequence of aligning by accelerometer 
leveling and gyrocompass ing. The attitude, velocity and position errors 
that result from these earth rate sensing errors may be computed with 
Equations A-l, A-4, A-5, A-ll, A-12, and A-13. 


A2 . 6 Accelerometer Bias. Level Accelerometer 

If the orientation of the IMU is not changed, an uncompensated level 
accelerometer bias will cause the following attitude, velocity and 
position errors: 

* r[’ ■ 005 v] (A - 18) 

Ep = - cos «> 0 tl (A-20) 

w o 
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where: Eg = level accelerometer bias, g 

A Y accelerometer bias error will result in a Z axis attitude error and 
Y axis velocity and position errors. Conversely, a Z accelerometer bias 
error will result in a Y axis attitude error and Z axis velocity and 
position errors. 

The attitude, velocity and position errors all contain periodic 
terms with period T . In addition, the attitude and position errors 
contain constant terms. 

The error sensitivity of the constant term in Equation A-20 ex- 
pressed in units of n.mi./yg is as follows: 



1 _ 


R 

g 


( 2 . 09x1 0 7 ) ( 1 0~ 6 ) 

6,000 


= .0035 H'- -- (A-21) 

yg 

where: E,-' = level accelerometer bias, wg 

o 

The attitude, velocity and position errors for an accelerometer 
bias error of 10 yg are plotted in Figures 3-7, 3-8 and 3-9. 

A2.7 Level Alignment Error 

An error in level alignment will cause a gravity computation error. 
The equivalent level accelerometer bias error is as follows: 

E b = E 6 g (fl-22) 


where: Eg = level accelerometer bias error, g 

Eg = level alignment error, rad. 

A Y axis alignment error will result in a Z accelerometer bias error. 
Conversely, a Z axis alignment error will result in a Y accelerometer bias 
error. 
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The accelerometer bias errors will cause attitude, velocity and 
position errors. These errors may be computed with Equations A-18, A-19 
and A-20. 

The error sensitivity of the constant position error term 
(Equation A-20) expressed in units of n. mi. /arc seconds is as follows: 

V _ (2.09x10 7 )(4.85x10~ 6 ) 

Eg 7- 6,000 

= .017 (A-23) 

arc sec 

where: Eg' = level alignment error, sec 

In the RSP , the reference coordinate system is level aligned by 
means of the accelerometers. Thus, a level accelerometer bias error will 
cause a level alignment error. This alignment error, in turn, will 
cause a gravity computation error that will offset the accelerometer bias 
error. Once navigation is initiated, this error offset will continue until 
the orientation of the IMU is changed. IMU rotation about a level axis 
will cause one of the two offsetting errors to change and thus cause the 
total error to increase. This offsetting error relationship is similar 
to that for azimuth alignment discussed in Section A2.3. 

In addition to the foregoing errors, a level alignment error will 
cause an earth rate transformation error. The error component in each of 
the level reference axes is as follows: 

E = E,- ft sin x (A-24) 

r o 

The corresponding errors in the vertical axis are as follows: 

E = E c ft cos X sin ^ (A-25) 

r o 

The earth rate transformation error will cause attitude, velocity, and 
position errors. The errors caused by the level axis components may be 
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computed with Equations A-l, A-4, and A-5. The errors caused by the 
vertical axis components may be computed with Equations A-ll, A-l 2, and 
A-l 3. 

The error sensitivity of the constant velocity error term in 
Equation A-4 expressed in units of (n.mi ./hr)/arc seconds is as follows: 
(x is assumed to be 35 deg). 

E ' 

= R q sin 35° 
o 

_ ( 2. 09x1 0 7 ) (4. 85x1 0~ 6 )( 3,600) (1 5) ( 4.85x1 0~ 6 )( .574) 

6,000 


= .0025 (A-26) 

arc sec 

The error sensitivity of the linear velocity error term in 
Equation A-l 2 evaluated at t = 45 min and expressed in units of 
(N-mi ./hr)/arc seconds is as follows: (x is assumed to be 35 deg and <p 

is assumed to be 90 deg). 

R q 2 (cos 35°) Z (sin 90°) 2 (45)(60) 

( 2 . 09x1 0 7 ) ( 4 . 85x1 0~ 6 ) C3 ,600 ) ( 1 5 ) 2 ( 4 . 85x1 0~ 6 ) 2 ( . 81 9 ) 2 ( 45 ) ( 60 ) 

6,000 

= .0006 n . - mi 7 h . r - (A-27) 

arc sec 

A2.8 Accelerometer Scale Factory Level Accelerometer 

A level accelerometer scale factor error will cause attitude, 
velocity and position errors that are a function of IMU velocity time 
history. For the special case of an acceleration impulse at the beginning 
of the navigation period and constant velocity (magnitude and direction) 
thereafter, the attitude, velocity and position errors are as follows: 

EyV sin to t 

<t> io R 

o 
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E 


(A-28) 



E y = E ? V cos u» 0 t 


(A-29) 


E = — sin o> t (A-30) 

p Uq o 

where: E^ = accelerometer scale factor error, dimensionless 

V = velocity increment, ft/sec 

In the forthcoming tests, the navigation errors caused by acceler- 
ometer scale factor errors will be much less than those caused by other 
error sources. 

A2.9 Accelerometer Misalignment, 

Level Accelerometers 

Misalignment of the level accelerometers with respect to the level 
gyros in the level plane will cause velocity errors that are a function 
of velocity time history. These errors are similar to the accelerometer 
scale factor errors discussed in Section A2.8. 

Misalignment of the level accelerometers with respect to the level 
gyros in the vergical plane will cause earth rate transformation errors. 

The effects of these errors are discussed in Section A2.7. 

A2.10 Vertical Velocity and Position Errors 

Vertical velocity and position errors will be limited by altitude 

damping in the RSP. It is assumed that the gain constants a and a will 

r v 

be set to -1 and 0, respectively. With these values, the vertical position 
error will be limited to the difference between the actual altitude and 
the altitude input, and the vertical velocity error will be small. 

A2.ll Navigation Errors of a Dodecahedron 

The navigation errors of a non-redundant strapdown IMU were dis- 
cussed in Sections A2.2 through A2.10. This discussion will now be 
extended to a dodecahedron. 

For the attitude sensors (SARs) of the BB DDH, 

m = Ab + e (A-31 ) 
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where: 


where: 


where: 


where: 


m = sensor output (6 vector), rad/sec 
A = transformation matrix (6x3 matrix), dimensionless 
= sensed quantity (3 vector), rad/sec 
£ = sensor error (6 vector), rad/sec 

b = B m (A-32) 

= least squares estimate of the quantity 
(3 vector), rad/sec 

B = (A T A) _1 A T 
b = b - b 

I) = error in the estimate of the sensed 
quanity (3 vector), rad/sec 

b = B £ 

The sensor error vector may be written as follows: 

e = C - M_W-,a + M T W,a + C,W,b 
— -o — s— l— —I Z— —l—o— 

+ C 2 W 4 b + ^Wgb (A-36) 

bias error (6 vector), rad/sec 

spin axis mass unbalance (6x6 diagonal 
matrix), rad/sec/g 

input axis mass unbalance (6x6 diagonal 
matrix), rad/sec/g 

sensor misalignment components (6x6 
diagonal matrices,), rad 

transformation matrices (6x3 matrices), 
dimensionless 

thrust acceleration (3 vector), g 
angular velocity (3 vector), rad/sec 


So 

H. 


fli 


-r -2* -3 


W r W 2 , Mg, = 

W 4 ’ ^5 


a_ = 
b = 


(A-33) 

(A-34) 

(A-35) 
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Assume that the level sensed acceleration components and the level 
angular velocity components are negligible. 



Equation A-36 may be rewritten as follows: 


* = ^ - M^a x + M I W'a x + 


where: W , Wp, W~ = transformation vectors (6x1 vectors), 

dimensionless 


(A-37) 


(A-38) 


' s ‘ 


C sin 3 a 


‘ C ' 

s 


C sin Bg 


C 

c 


S sin B r 


S 

c 

— 2 = " 

L 

S sin Bp 

w 1 = 

-3 

S 

0 


sin b e 


1 

0 


sin Bp 


1 


where: S = sin a 


C = cos a 
a = .55 rad 

6^ = gimbal angle of the ith SAR (the SAR gimbal angle 

is measured from the position in which the OA 
is in the plane of the IA and the DDH X axis 
and is in the positive X direction), deg 


(A-39) 
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The sensor misalignment matrix is 




HA 

0 

0 

0 

0 

0 


0 

: 1B 

0 

0 

0 

0 


0 

0 

: ic 

0 

0 

0 


0 

0 

0 

HD 

0 


0 

0 

0 

0 

HE 

0 


0 

0 

0 

■ 0 
0 
c 


(A-40) 


IF 


where: C 


li 

! 0i 

■n 


'oi + 'li s1n S i 

external misalignment of the ith SAR, rad 
internal misalignment of the ith SAR, rad 


With no failures, the transformation matrix in Equation A-31 is 


A = 


S 

S 

c 

c 

0 

0 


The B matrix (Equation A-32) is thus 


B = (AV 1 A T 


0 

0 

-S 

S 

C 

C 


C 

-C 

0 

0 

S 

-S 


(A-41) 


S S C 

0 0 -S 

C- -C 0 


C 0 

s c 
0 s 


0 

c 

-s 


(A-42) 


Assume now that the uncompensated SAR biases are statistically 

2 

independent with zero mean and variance equal to Og . Thus 


- a r 


= Or 


= O r 


A 


(A-43) 


A- 18 



Considering only the bias term in Equation A-38, the sensor error co- 
variance matrix is 


T 

<£ £ > 


h 


where: 1^ = 6x6 identity matrix. 

The covariance matrix of the error in the estimate of the sensed quantity 
is 


<b b T > = <B ^ £ T B T > 

= <(A T A) _1 A T £ £ T a(a t a) -1 > 


or 


= (a t a) _ 1 Og i 3 


°B I 


’ a b 

X 


.707 

a b 

= °B 

.707 

y 


,° b 2 . 


.707 


2 

where: a. = variance of the ith component of the error 

D i in the estimate of the sensed quantity. - 


Similarly, for the spin axis mass unbalance term. 


T 

< E £ > 


.. 2 2 
" a s a x 


s 2 0 0 0 

0 s 2 0 0 

0 0 c 2 0 

000c 2 
0 0 0 0 


0 0 
0 0 
0 0 
0 0 
0 0 


0 0 0 


0 0 


0 


(A-44) 
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2 

where: o s 


= spin axis mass unbalance variance 


2 2 
a a 

; ;T _ s x 
<b b > = — 3 — 


2(S 4 + C 4 ) 
0 


2S 2 C 2 


0 

0 


and 


<e e > = a 


where: a 


I 



l 

0 

0 


2S 

2 C 2 

’ °b 

X 



.548 




a b 

= a 

a 

.316 




y 

S 

X 





a b 



.316 




z 



- 




is mass 

unbalance 

term with 



6 A = 0 

B 

- = B 

P = 90 deg 




r 2 

(T 

0 

0 0 


0 

0 


0 

c 2 

0 0 


0 

0 


0 

0 

S 2 0 


0 

0 

2 2 







I a x 

0 

0 

0 S 2 


0 

0 


0 

0 

0 0 


1 

0 


_ 0 

0 

0 0 


0 

1 

s mass 

unbalance 
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<b b > = 


2 

4S 2 C 2 


0 

0 

a x 

4 

0 

2(S 4 

+ c 4 ) 

0 


0 


0 

2(S 4 + C 4 )_ 

1 

Q 

cr 

X 



.447 


°b 

y 

= °I 

a x 

.548 


. V 



.548 
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For the SAR misalignment term, and ignoring the internal mis- 
alignments 


’ °b 


-P* 
1 

X 



°b 

= a b 

.548 

y 

c X 


“i 

N 

-Q 

O 

i 


.548 
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A3. 0 FAILURE DETECTION, DIAGNOSIS AND CORRECTION 

A failure of one or more sensors will cause system performance to 
be degraded. This degradation will be discussed in terms of the following 
effects: 

a. Undetected errors 

b. Transient errors 

c. Degradation resulting from the deletion of sensors. 

Undetected errors are failure induced errors that do not exceed the 
detection threshold. Transient errors are errors that develop in the 
interval from the occurrence of the failure to its detection and correction. 
The effects of transient errors may remain after the failure is cleared. 

The degradation resulting from the deletion of sensors refers to the 
increased error resulting from the deletion of sensors in the estimation 
equations. 

Undetected errors and transient errors are both a function of the 
filter gain and time constant values in the FDDC subroutine of the RSP. 
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The effects of undetected errors may be evaluated with the equations in 
Section A2.0 and the FDDC text matrix (C matrix) in Reference 3. The 
effects of transient errors may be evaluated with the RSP in the Simulation 
Mode. (The effects of undetected errors may also be evaluated in the 
Simulation Mode). 

The performance degradation resulting from the deletion of sensors 
may be evaluated with the equations in Section A2.ll. The effect of one 
or more sensor failures may be determined by setting the appropriate row 
(or rows) in the A matrix (Equation A-41) equal to zero. 

A4.0 ESTIMATED PERFORMANCE 

Navigation accuracy in the forthcoming tests may be estimated from 
the equations in Section A2.0 and the estimated performance characteristics 
of the BB DDH. The latter information was supplied by MSFC and is sum- 
marized in Table A-l . 

The navigation error estimates will be helpful in determining 
whether or not particular RSP changes are required. These estimates will 
also be helpful in designing the laboratory and van tests and in evaluating 
the results of these tests. 

Five cases will be considered: 

a. DDH not rotated during the navigation period, 
no prenavigation gyro drift update, gyro com- 
pensation not a function of SAR gimbal angle. 

b. DDH not rotated during the navigation period, 
gyro drift updated prior to navigation, gyro 
compensation not a function of SAR gimbal angle. 

c. DDH rotated about the vertical axis at the 
beginning of the navigation period, no prenavi- 
gation gyro drift update, gyro compensation not 
a function of SAR gimbal angle. 

d. DDH rotated about the vertical axis at the 
beginning of the navigation period, gyro drift 
updated prior to navigation, gyro compensation 
not a function of SAR gimbal angle. 
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Table A-l 

DDH Performance Characteristics 


Error Source 


Gyro Bias 


Day to day stability 
Single run stability 


Gyro Spin Axis Mass Unbalance 


Gyro Misalignment 
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Day to day stability 

.05 deg/hr/g 

Single run stability 

.01 deg/hr/g 

Gyro Input Axis Mass Unbalance 

Magnitude 

.10 deg/hr/g 

Day to day stability 

.05 deg/hr/g 

Single run stability 

.01 deg/hr/g 

Gyro Output Axis Unbalance 

Magnitude 
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Day to day stability 

Single run stability 
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e. DDh rotated about the vertical axis at the beginning 
of the navigation period, gyro drift updated 
prior to navigation, gyro mass unbalance and mis- 
alignment compensation both a function of SAR 
gimbal angle. 

Only the nonperiodic error components will be considered in the estimates. 
Dynamic errors and the performance degradation caused by simulated 
failures will not be considered. 

CASE A 

DDH not rotated during the navigation period, no prenavigation 
gyro drift update, gyro compensation not a function of SAR 
gimbal angle. 

1. Azimuth alignment 

Error contribution negligible. 

2. Level alignment 

(0025) (10x1 0 -6 ) (2. 06x1 0 5 ) @ ( .0006) (100x1 O' 6 ) 

(2.06xl0 5 ) « .05 n.mi./hr 

3. North axis attitude drift 

60.8[( . 1 0) ( .707) @ ( .05) ( .316) Q ( . 02 ) ( .548) J = 4.45 n.mi./hr 

4. East axis attitude drift 

60.8[( .02) ( .707) © ( .01 ){ .316) © ( .01 )( .548)3 = .94 n.mi ./hr 

5. Vertical axis attitude drift 

9.7[(.10)( .707) © ( . 05 ) ( . 548 ) © (.02)(.447)] = .74 n.mi ./hr 

6. Gyro misalignment 

(.0025) (15) (.548) © ( .0006) (15) ( .447) » .02 n.mi./hr 

7. Accelerometer bias 

Error contribution negligible. 

8. Accelerometer scale factor 
No nonperiodic error term. 
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9. Accelerometer misalignment 

( . 0025) (1 5) Q (.0006) (15) w .04 n.mi./hr 

10. Total error 

.05 0 4.45 0 .94 0 .74 0 .02 0 .04 « 4.6 n.mi./hr 
CASE B 

DDH not rotated during the navigation period, gyro drift updated 
prior to navigation, gyro compensation not a function of SAR 
gimbal angle. 

1. Azimuth alignment 

Error contribution negligible. 

2. Level alignment 

Error contribution negligible. 

3. North axis attitude drift 

60. 8[ ( .02) ( .707 ) 0 (.01 )( .316) 0 (.01 )( .548)] = .94 n.mi./hr 

4. East axis attitude drift 

60.8[ ( .02) ( .707 ) 0 ( .01 ){ .316) 0 ( .01 )( .548)] = .94 n.mi./hr 

5. Vertical axis attitude drift 

9 . 7 [ ( .02) ( . 707 ) 0 ( .01 ) ( .548) 0 ( .01 )( .447)] = .15 n.mi ./hr 

6. Gyro misal ignment 

Error contribution negligible. 

7. Accelerometer bias 

Error contribution negligible. 

8. Accelerometer scale factor 
No nonperiodic component. 

9. Accelerometer misalignment 
Error contribution negligible. 

10. Total error 

.94 (0 .94 (0 .15 *1.3 n.mi ./hr 


A-25 



CASE C 


DDH rotated about the vertical axis at the beginning of the 
navigation period, no prenavigation gyro drift update, gyro com- 
pensation not a function of SAR gimbal angle. 

1. Azimuth alignment 

60 . 8[ ( . 1 0 ) ( . 707 ) © ( .05) ( .316) © ( .02)( .548)] = 4.45 n.mi./hr 

2. Level alignment 

( .0025) (100x1 0 -6 ) (2 .06x1 0^) © ( .0006) (100x1 0 -6 ) 

(2.06x10^) « .05 n.mi./hr 

3. North axis attitude drift 

60.8[( . 1 0 ) ( .707) © ( . 05) ( . 31 6 ) © ( -02) ( .548) © 

( .10) ( .276)] = 4.78 n.mi./hr 

4. East axis attitude drift 

60.8[( . 10) ( . 707 ) © ( .05) ( .316) © (.02)(.548) © 

( .10) ( .724)] = 6.25 n.mi./hr ^ 

5. Vertical axis attitude drift 

9.7[( .10)( .707) © ( . 05 ) ( .548) © ( . 02 ) ( .447) © 

( . 1 0) ( . 90)3 = 1.14 n.mi./hr 

6. Gyro misalignment (assume 360 deg rotation) 

.106[(30)( .276) © (30) ( .724)] « 2.46 n.mi. 

. 0036[ (15) ( .447 ) © (30) ( . 90)] © .0025[(1 .73)(15)( .548) 

© (30) ( . 276) © (30)( .724)] © .0006[(15)( .548) © (15) ( .447) 
© ( .707) (30) ( .276) © ( . 707) (30) ( . 724)] « .12 n.mi ./hr 

7. Accelerometer bias 

Error contribution negligible. 

8. Accelerometer scale factor 
No nonperiodic component. 

9. Accelerometer misal ignment 

. 0025[ 1 5 © 15] © .0006[1 5] = .05 n.mi./hr 
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10. Total error 


2.5 n. mi 

4.45 © .05 © 4.78 © 6.25 © 1.14© .12 © .05=9.1 n/ni./hr 
CASE 0 

DDH rotated about the vertical axis at the beginning of the nav- 
igation period, gyro drift updated prior to navigation, gyro compen- 
sation not a function of SAR gimbal angle. 

1. Azimuth alignment (assume 3 min misalignment) 

60.8)[(15) ( .81 9) (3) (60) (4.85x1 0 -6 ) © (.02)(.707) © (.01)(.316) 
© ( .01 ) ( .548)] = 1.14 n.mi ./hr 

2. Level alignment 

Error contribution negligible. 

3. North axis attitude drift 

60.8[( .02) ( .707) © ( .01 ) ( .316) © (.01)(.548) © 

( . 10) ( .276)] = 1.92 n.mi. /hr 

4. East axis attitude drift 

60 . 8[ ( 1 5 ) ( .819) (3) (60) (4. 85x1 0" 6 ) © ( .02) ( .707) © 

( .01 ) ( .316) © ( .01 )( .548) © ( .10)( .724)] = 4.54 n.mi. /hr 

5. Vertical axis attitude drift 

9.7[( .02) ( .707 ) © ( .01 )( .548) © (.01)(.447) © 

(.10)(.90)] = .88 n.mi. /hr 

6. Gyro misalignment (assume 360 deg rotation) 

. 1 06[ ( 30 ) ( . 276 ) © (30) ( . 724) ] « 2.46 n.mi 

.0036[(15)( .447) © (30)(.90)] © .0025[(1 .41 ) (15) ( .548) © 
(30)(.276) © (30) ( .724)] © . 0006[ ( 1 5) ( . 548) © ( .707) (30) 

(.276) (T) ( . 707 ) (30) ( .724)3 « .12 n.mi. /hr 

7. Accelerometer bias 

Error contribution negligible. 
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8. Accelerometer scale factor 
No non-periodic component. 

9. Accelerometer misalignment 
(.0006) (15) = .01 n.mi./hr 

10. Total error 
2.5 n.mi . 

1.14 © 1.92 Q 4.54 © .88 © .12 © .01 « 5.1 n.mi./hr 
CASE E 

DDH rotated about the vertical axis at the beginning of the 
navigation period, gyro drift updated prior to navigation, gyro mass 
unbalance and misalignment compensation both a function of SAR 
gimbal angle. 

1. Azimuth alignment (assume 3 min misalignment) 

60. 8[ ( 1 5) ( .81 9) (3) (60) (4.85x1 0 -6 ) © (.02)(.707) © 

( .01 ) ( .316) © ( .01 )( .548)] = 1.14 n.mi ./hr 

2. Level alignment 

Error contribution negligible. 

3. North axis attitude drift 

60.8[( .02) ( .707 ) © ( .01 )( .316) © ( .01 )( .548)] = .93 n.mi ./hr 

4. East axis attitude drift 

60.8[(15) ( .819) (3) (60) (4.85xl0~ 6 ) © (.02) (.707) © 

( .01 )( .316) © ( .01 )( .548)] =1.14 n.mi./hr 

5. Vertical axis attitude drift 

9.7[( .02)( .707) © ( .01 )( .548) © ( . 01 )( .447)] = .15 n.mi ./hr 

6. Gyro misalignment 

. 10611(15) ( .276) © (1 5) ( .724)] 1.23 n.mi. 

. 0036[ ( 1 5 ) ( .447)] © .0025[(1 .41 )(15)( .548)] © 

. 0006 C ( 1 5 ) ( . 548 ) ] « .04 n.mi./hr 
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7. Accelerometer bias 


Error contribution negligible. 

8. Accelerometer scale factor 
No nonperiodic component. 

9. Accelerometer misalignment 
( . 0006 ) ( 1 5 ) = .01 n.mi./hr 

10. Total error 
1 .2 n.mi . 

1.14 @ .94 0 1.14 Q .15 0 .04 0 .01 ~ 1.9 n.mi ./hr 
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APPENDIX B 

Bl. NOMINAL NAVIGATION RUN 


The following tables contain the output word list and 6-sec printouts 
from a test case dated 25 May 1969. 
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The following tables contain the output word list and 3-sec printouts 
from a test case dated 28 May 1969. 
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Appendix C 

FDDC FILTER CONSTANTS 


The method of deriving initial values of the FDDC prefilter gain and 
time constant is given here. 

Cl .0 FILTER DESCRIPTION 

I 

V is a vector of 15 test signals each of the form: 

V^ n = CXI + 0X2 - SY1 + 0Y2 - CZ1 + SZ2 

where C and S are the sin and cosine of a, (.55 radian) 

and XI, X2, Y1 , Y2, Z1 , Z2 are the sensor outputs. 

Before comparison to the threshold value (unity), each of the 15 gyro 
test signals is passed through a first order low pass filter of the form: 


Fi = 


T f 


Kf »' 
~S+T V in 


(Cl) 


If V. is a step function 

Fi = V.' n Kf (l-e" t/T ) (C2) 

= si net we Kf (l-e" t/,T ) 

where Fi > 1 constitutes a failure state 
and we is the failed gyro error (bias) 

Kf is the filter gain 
t is the filter time constant 

and a is the DDH alignment angle, .55 radian. 

★ 

or COSa 
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C2.0 FILTER THRESHOLD 


The steady state filter threshold is 1/Kf, computed as follows: 
For t>>T, equation (C2) becomes: 


Fi - »,„ Kf 


Since Fi = 1 constitutes a failure, the threshold (for V^) is: 


TH V = -j^r (radian/second) 


(C3) 


The failure threshold for the output of a single gyro is 1/Kf sina, or: 


TH = 


1.9 


g Kf 


(C4) 


C3.0 QUANTIZATION VERSUS GAIN AND TIME CONSTANT 

For intervals short compared to the time constant 
For t<<T, equation (C2) becomes: 


Fi = V 


in " (f) 


(C5) 


At failure Fi = 1 , therefore, the maximum rate is: 


V in (max) ■ r 


Kf 


(C6) 


or, expressed in terms of error angle 


’(max) = V in (max)t = W - ™ (t) 


(C7 ) 


for periods short compared to the time constant. 
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The maximum quantization as a function of filter parameters is, 
therefore , 


<U * W < C8 > 

Or, given the necessary filter threshold and the quantization, the 

minimum permissible time constant, t , is: 

m 

■ <5 Kf - ?H < c9 > 

C4.0 COMPUTING TIME TO FAILURE 

From (C2), failure is indicated when 

Vi Kf (1 -e -t/x ) = 1 (CIO) 

Solving for t 



C5.0 COMPUTATION OF FILTER CONSTANTS 

The minimum failure threshold that will not trigger on normal (3a) 
instrument behavior plus noise is desired. 

First an estimate must be made of the expected total uncertainty 
in each gyro error (fixed drift @ g sensitive @ misalignment error X 
constant rate (+) dynamic errors) under actual test conditions, after 
instrument compensation. Then the expected variance in the test signal 
V i should be computed by adding the variances of four gyros (but 
neglecting quantization and high frequency noise), and finally compute 
aV i, the (expected) standard deviation of the test signals. 
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C6.0 SELECTION OF THRESHOLD (GAIN) 


Some margin must be allowed between expected instrument .behavior and 
the failure threshold. The necessary margin is a function of the noise 
present and the degree of filtering. For a first assumption the initial 
value of the failure threshold will be set at 4 aV^. 

From (C3) 


Kf = PT (second/radian) (C12) 

I 

In estimating aV.. it is seen that if constant high rate turns 
(10 degrees/second) are allowed, and a total misalignment uncertainty of 
1 arc minute, 3a is assumed, 4aV.. is on the order of 
sin (1 arc minute) X 36,000 = 12 degrees/hour. Setting the threshold at 
this level would prohibit detection of small but intolerable bias shifts. 

As an alternative, if we neglect the misalignment error in 
estimating 4aV i , a preliminary estimate is 0.2 degrees/hour or 

TH = .2 degrees/hour — 10" 6 radian/second (Cl 3 ) 

and 



second/radian 


(C14) 


C7.0 SELECTIONS OF TIME CONSTANT x 
From (C9) 


= Q Kf (CIS) 

For the BB DDH the gyro quantization, A 0 , = 2.4 x 10” 5 radian and a 

I c 

worst case estimate of the quantization on V.. is 7 x 10 radian 

.*• = (7 x IQ* 5 ) Kf = 70 seconds (Cl 6 ) 
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where x m is the minimum time constant to avoid triggering the failure 
detection logic on quantization noise. 

For t five to ten times larger than t , quantization noise is 
effectively suppressed at the expense of transient response. The initial 
value of t will be set at 500 seconds, the time used in digital simula- 
tions. 

C8.0 EVALUATION OF SELECTED CONSTANTS 
Using the values of 

Kf = 10^ second/radian 


and 


t = 500 seconds 


from (C7) 


max 


Kf 


^jr radian - 103 arc second 
10 6 


(Cl 7) 


Again assuming a sensor misalignment angle y, of 1 arc minute, the 
maximum permissible vehicle turn (in a time short relative to the time 
constant, e.g., < 150 seconds) would be (approximately) 


0 . 


max 


sina siny 


-4 

5.10 radian 
.00015 


= 190 degrees 


(Cl 8 ) 


The required threshold to allow a 360 degree turn would be 
.38 degrees/hour which is large relative to permissible gyro drift errors. 

The selected values of Kf = 10^ and t = 500 are, therefore, 
considered adequate for initial values, subject to the limitations dis- 
cussed; turns will be limited probably to 180 degrees, by the misalign- 
ment uncertainty. 
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C9.0 TEST CASE RESULTS 

Table C-I is a summary of results of a filter simulation using the 
recommended values of: 

t = 500 seconds 
Kf = 10^ second/radian 
Threshold = .2 degrees/hour 

and quantization of the test signal (summation of 4 gyro outputs) is 
13 arc seconds. 

Table C-I. Filter Simulation Results 


Gyro Test/, x 
Signal Error ' ' 

°/hr (1) 

Predicted 

Time of^ 
Failure-sec 



Actual' ' Time 
of 

Failure-sec 

Accumulated 
Error Angle 
Arc seconds 

.2 

oo 


oo 

.25 

804 

725 

182 

.4 

346 

329 

143 

.8 

144 

162 

117 

1.6 

67 

66 

117 

3.2 

32 

32 

104 

6.4 

16 

17 

104 

10 

10.1 

9.3 

104 

12.8 

7.8 

8.5 

104 

25.6 

3.9 

4.4 

104 

250 

.4 

.44 

104 


^Corresponds to sina (gyro error) or .53 (gyro error) 

(?) 

' 'Differs from predicted time due to random quantization error 


Figure C-l is a plot of failure detection time and accumulated 

I 

error angle versus V ../threshold , the ratio of the test signal error to 
the threshold. 
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Cl 0.0 COMPUTATION OF DIGITAL FILTER CONSTANTS 


The digital form of the prefilter is 


F (ir Kf l F (1-l) + Kf 2 V in 


where V.. is the test signal, (sina we) 


Kf i = e 


-T/t 


where T is the sample interval, .04 seconds 


and 


Kf 2 = Kf (1 -Kf 1 ) 


For t = 500 seconds 


Kf -j = e "- 04/50 ° = .99992 

g 


Substituting (Cl 4 ) and (C22) in (C21) 


Kf^ = (1 -.99992) 10° = 80 seconds/radian 


C11.0 SUMMARY 

The initial values of the FDDC prefilter are, therefore: 


Kf ] = .99992 
9 

Kfg = 80 seconds/radian 


For a filter threshold of .2 degrees/hour (on V.) and a time 
500 seconds. 
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( Cl 9 ) 


(C20) 


(C21 ) 


(C22 ) 


(C23) 


constant of 



